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ABSTRACT 

The optical spectrum and light curve of EC13471-1258 shows that it is an eclipsing 
binary with an orbital period of 3^ 37 m comprising a DA white dwarf and a dMe 
dwarf. Total eclipses of the white dwarf are observed lasting 14 min, with the partial 
' phases lasting 54 s. On one occasion, two pre-eclipse dips were seen. Timings of the 

, eclipses over ten years show jitter of up to 12 s. Flares from the M dwarf are regularly 

^5 observed. The M dwarf also shows a large amplitude ellipsoidal modulation in the V 

band light curve. The component stars emit almost equal amounts of light at 5500 A. 
HST STIS spectra show strong Lyman alpha absorption with weak metal lines of 
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Q ■ C I, II and Si II superimposed. Model atmosphere analysis yielded an effective tem- 

£j I perature of 14220 ± 300 K and log g of 8.34 ± 0.20 for the white dwarf with these 

errors strongly correlated. Its metal abundance is l/30th solar with an uncertainty of 
0.5 dex, and it is rapidly rotating with Vi sin i = 400 ± 100 km s _1 . The white dwarf 
also shows radial velocity variations with a semi-amplitude of 138 ± 10 km s _1 . The 
gravitational redshift of the white dwarf was measured: 62 km s . 

From optical spectroscopy the spectral type of the M dwarf was found to be M3.5- 
M4, its temperature 3100 ± 75 K, its rotational velocity 140 ± 10 km s _1 , its radial 
velocity semi-amplitude 266 ± 5 km s _1 , its mean V — I colour 2.86 and its absolute 
V magnitude 11.82. Intriguingly, its metal abundance is normal solar. 

The Ha emission line shows at least two distinct components, one of which is 
uniformly distributed around the centre of mass of the M dwarf and provided the esti- 
mate of the rotational velocity of the M dwarf. The other arises from the other side of 
the binary centre of mass, well within the white dwarf Roche lobe. This behaviour is 
confirmed by Doppler tomography which shows the presence of two distinct velocity 
components within the primary Roche lobe. The interpretation of these features is un- 
certain. Variations in strength of the components with binary phase can be attributed 
to optical thickness in the Balmer lines. Similar behaviour is seen in the observations 
of the other Balmer emission lines, although with poorer signal-to-noise. Flares in Ha 
were observed and are consistent with arising from the vicinity of the M dwarf. 

Dynamical solutions for the binary are discussed and yield an inclination of 75.5 ± 
2.0°, a white dwarf mass and radius of 0.78 ± 0.04 M and 0.011 ± 0.01 R , and an 
M dwarf mass and radius of 0.43 ± 0.04 M and 0.42 ± 0.02 R Q . These parameters 
are consistent with the Wood (1995) mass-radius relation for white dwarfs and the 
Clemens et al. (1998) mass-radius relation for M dwarfs; we argue that the M dwarf 
just fills its Roche lobe. The radius of the white dwarf and the model fit imply a 
distance of 48 ± 5 pc and an absolute V magnitude of 11.74. 

The rapid rotation of the white dwarf strongly suggests that the system has un- 
dergone mass transfer in the past, and implies that it is a hibernating cataclysmic 
variable. The M dwarf shows the properties expected of secondaries in cataclysmic 
variables: chromospheric activity and angular momentum loss. 
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1 INTRODUCTION AND SYSTEM 
OVERVIEW 

Stellar evolution is one of the most secure theories in mod- 
ern astrophysics; masses and radii of stars are fundamental 
parameters of the stars to which this theory applies. When 
stars are found in interacting binaries, additional complex- 
ities are introduced; results from single star evolution are 
often applied for want of any alternative approach. Systems 
which offer the prospect of testing these prejudices are there- 
fore worthy of some scrutiny. 

Prominent among the interacting binaries are cata- 
clysmic variables: short period binaries comprising a Roche 
lobe filling late K or M dwarf and a white dwarf. If the white 
dwarf's magnetic field is small, the mass that is transferred 
from the cool star to the white dwarf forms an accretion 
disc. The energy generated by the accretion almost always 
far exceeds the luminosity of the component stars, afford- 
ing a valuable laboratory for studying accretion processes. 
However, masses, radii, temperatures and other fundamental 
properties of the component stars are usually very difficult 
to determine because of the glare of the accretion light. 

The progenitors of cataclysmic variables are detached, 
short period DA+dM binaries which have emerged from 
common envelope evolution. This process is thought to give 
rise to all very short period binaries and although there has 
been a lot of theoretical work in this area, comprehensive 
observational tests of the theory are required. 

In this paper we report an observational study of a DA 
+ dMe eclipsing binary with an orbital period of 3 h 37 m . 
As we will show, this system is almost a cataclysmic vari- 
able: there is a little mass transfer taking place, but the 
component stars can nevertheless be studied with relative 
ease, their fundamental properties derived and compared 
with theory. 

EC13471-1258 was discovered as part of 
the Edinburgh-Cape blue object survey (Stobie et al. 1997, 
Kilkenny et al. 1997). Eclipses were detected serendipitously 
by one of us (DOD) during measurement of the U BV colours 
of the object: the observer could not identify why the count 
rate from the star declined suddenly, despite the apparently 
clear weather and normal behaviour of the instrument. The 
subsequent "unassisted" recovery of the count rate heralded 
the realization that a rapid eclipse had been observed and a 
program of high speed photometry was initiated to study the 
eclipses. These data will be considered in detail in the next 
section. Fig. 1 (upper) shows a Johnson V band light curve 
at 30 s time resolution. The binary period, 0.151 d (3.62 hr), 
can be easily estimated from the horizontal axis. Also obvi- 
ous is a modulation at half the orbital period resulting from 
the observer viewing at different aspect the geometrically 
distorted uneclipsed star (an "ellipsoidal modulation"). Fig. 
1 (lower) shows one of the better "white light" eclipses at 
1 s time resolution. The eclipse is 1.5 mag deep and total. 
The time from first to last contact is 944 s, the duration of 
totality is 835 s and ingress (or egress) take place in 54 s. 



Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract NAS 5-26555. 
These observations are associated with program 7744. 
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Figure 1. (Upper) V band light curve at 30 s time resolution. 
The ordinate is extinction corrected counts s _1 in the natural 
unfiltercd, photometric system. (Lower) White light eclipse light 
curve at 1 s time resolution. The ordinate is extinction corrected 
counts s _1 . 
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Figure 2. Spectrum of EC13471-1258 on 1993 Apr 15. The ab- 
scissa is in A and the ordinate is in erg s — 1 cm~ 2 A -1 . The ac- 
curacy of the flux distribution shortward of 3800 A and longward 
of 6800 A is uncertain. 



The eclipses are slightly deeper in white light than in the V 
band. 

The absence of a secondary eclipse indicates that the 
cool secondary is much larger than the hot primary; this is 
confirmed by the optical spectrum of EC13471-1258 which is 
shown in Fig. 2 at 7 A resolution. These data were obtained 
through a wide slit in order to obtain an accurate estimate 
of the flux distribution; we believe that the raised flux below 
3800 A and sagging flux above 6800 A are instrumental arti- 
facts, but that the rest of the spectrum is accurate. The flux 
distribution is quite flat: in the red region of the spectrum, 
the TiO molecular bands and other features typical of an M 
dwarf are visible, while the most prominent features in the 
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Table 1. High speed white light photometry and eclipse timings of EC13471-1258 



Run 


Date 


Start of Run 


Length 


Cycle 


Mid Ingress 


Mid Egress 


Name 




JD 


(hr) 




JDq 


JDq 






2440000+ 






2440000+ 


2440000+ 


S5465 


06 3 1992 


8687.553399 


2.1 








S5469 


07 3 1992 


8689.460013 


1.7 








S5470 


08 3 1992 


8689.570206 


1.9 








S5473 


10 3 1992 


8691.585272 


0.6 


13 


8691.595195 




ck0008 


02 4 1992 


8714.511519 


0.3 


165 




8714.520633 


ck009b 


02 4 1992 


8715.409417 


0.5 


171 


8715.414876 


8715.425178 


ckOOll 


03 4 1992 


8715.557234 


0.6 


172 


8715.565678 


8715.575927 


S5478 


09 4 1992 


8721.587396 


0.5 


212 


8721.595947 


8721.606231 


S5480 


11 4 1992 


8723.551101 


0.4 


225 


8723.555775 


8723.566070 


S5494 


27 5 1992 


8770.222061 


2.2 


535 


8770.290602 


8770.300917 


S5494a 


27 5 1992 


8770.315458 


2.0 








S5494b 


27 5 1992 


8770.398786 


1.4 








S5498 


28 5 1992 


8771.338662 


0.5 


542 


8771.345941 


8771.356243 


S5585 


24 2 1993 


9043.449210 


0.7 


2347 


9043.463426 


9043.473720 


S5587 


25 2 1993 


9043.607556 


0.5 


2348 


9043.614198 




S5591 


26 2 1993 


9044.509713 


0.6 


2354 


9044.518725 


9044.529037 


S5597 


27 2 1993 


9046.463340 


0.8 


2367 


9046.478580 


9046.488874 


ck29 


194 1993 


9097.423855 


0.6 


2705 


9097.434649 


9097.444946 


S5621 


21 6 1993 


9160.288209 


0.6 


3122 


9160.300567 


9160.310872 


S5678 


15 1 1994 


9367.581451 


0.6 


4497 


9367.592196 


9367.602487 


S5702 


10 3 1994 


9421.551509 


0.6 


4855 


9421.563384 


9421.573684 


S5704 


11 3 1994 


9422.605664 


0.6 


4862 




9422.628966 


ck70 


16 5 1994 


9489.395598 


2.3 


5305 


9489.404255 




S5754 


10 6 1994 


9514.406446 


0.9 


5471 


9514.430010 


9514.440315 


S5760 


12 6 1994 


9516.226818 


0.6 


5483 




9516.249406 


S5831 


03 3 1995 


9779.605904 


0.5 


7230 


9779.612423 


9779.622709 


S5837 


04 3 1995 


9781.411594 


0.3 


7242 


9781.421522 




S5840 


05 3 1995 


9782.469310 


0.5 


7249 


9782.476814 


9782.487129 


S5853 


24 5 1995 


9862.220185 


0.6 


7778 


9862.227543 


9862.237838 


S5858 


27 5 1995 


9865.392234 


0.4 


7799 




9865.403793 


ckl20 


31 5 1995 


9869.405236 


1.7 


7826 


9869.463916 


9869.474208 


ckl22 


01 6 1995 


9870.215706 


0.3 


7831 


9870.217693 


9870.227990 


dmk024 


28 1 1996 


10110.515827 


0.5 


9425 


10110.525155 


10110.535446 


S5865 


22 2 1996 


10135.542653 


0.5 


9591 


10135.550900 


10135.561204 


S5868 


23 2 1996 


10137.495595 


0.4 


9604 


1/11 o f" -| f\*-T A ~\ 

10137.510741 




S5871 


25 2 1996 


10138.560961 


0.4 


9611 


10138.566058 


10138.576355 


m0165c 


22 3 1996 


10164.548170 


0.4 








m0280b 


15 7 1996 


10280.237416 


1.3 


10551 


10280.278120 


10280.288420 


S5892 


14 2 1997 


10493.585765 


0.7 


11966 


10493.600012 


10493.610315 


dodOOl 


06 5 1997 


10575.305595 


0.5 


12508 


10575.310609 


10575.320907 


dod025 


19 8 1998 


11045.216518 


0.4 


15625 


11045.221857 


11045.232149 


ml254d 


17 3 1999 


11254.558177 


1.9 


17014 


11254.623987 


11254.634289 


runOOl 


12 5 1999 


11311.386676 


2.1 


17391 


11311.459600 


11311.469906 


dod052 


21 1 2002 


12295.600757 


0.4 


23919 


12295.604485 


12295.614793 



blue are the broad Balmer absorption lines of a white dwarf. 
In addition, narrow emission cores to the Balmer lines can 
be seen (especially at Ha and H/3) as well as the Call H 
and K lines in emission. 

Figs. 1 and 2 show that EC13471-1258 is a short pe- 
riod eclipsing binary comprising a white dwarf and an M 
dwarf. The plan of the paper is first to present all the ob- 



servational data and associated analysis. The last part will 
use the observational constraints to determine the binary 
system parameters and discuss the evolutionary status of 
the system. The approach will be similar to that applied to 
V471 Tau by O'Brien, Bond & Sion (2001). An earlier ob- 
servational study of EC13471-1258 has been published by 
Kawka et al. (2002). 
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2 HIGH SPEED WHITE LIGHT 
PHOTOMETRY 

Time series photometry was obtained in "white light" with 1 
s integrations using blue sensitive photomultiplier tubes on 
the SAAO 0.75 and 1.0 m telescopes. Data were obtained 
during nine observing seasons spanning ten years. An ob- 
serving log is given in Table 1. 

The photomultiplier data were reduced by subtracting 
the sky background and correcting for atmospheric extinc- 
tion using a mean extinction coefficient of 0.3 mag/airmass. 
Four different photomultiplier tubes were used; although all 
were blue sensitive, small differences in red cutoff may affect 
the depth of the eclipses. This possibility was kept in mind in 
the analysis presented below. Care was taken to ensure that 
the heliocentric Julian date was calculated as accurately as 
possible so that the eclipse times are precisely determined. 

We searched the longer light curves for evidence for 
rapid oscillations. No significant peaks in the Fourier am- 
plitude spectrum were found, with a limit of 1 mmag. 



2.1 The ephemeris 

The best quality light curves were selected and times of 
mid ingress and mid egress were measured by fitting three 
straight lines to the data around ingress (or egress) by least 
squares. One line was fitted to the out-of-eclipse data (with 
zero slope, of course). A second line with zero slope was 
fitted to the section in totality just after ingress (or before 
egress) . The third straight line was fitted to the ingress (or 
egress) section. The time of mid ingress (or egress) was taken 
to be the time when the third line crossed the average levels 
of the first two lines. These timings (for both ingress and 
egress are collected in Table 1). Run S5470 defines cycle 
zero but no timings are given. The reason for this is that 
the eclipse observed in S5470 was distorted by thin cirrus. 
Cloud interference usually explains the missing values as- 
sociated with some of the other runs in Table 1 (including 
those where only the ingress or egress timing is listed). 

By the third observing season it became apparent that 
the eclipses could not be explained by a linear ephemeris (i.e. 
a constant period). A quadratic ephemeris was thus fitted 
to the eclipse timings as more data were accumulated for a 
few years. However, in the last several years this quadratic 
ephemeris has proved to have little predictive power. Figs. 
3(a) and (b) show the timing residuals for egress and ingress, 
respectively, with respect to linear (Fig. 3a) and quadratic 
(Fig. 3b) ephemerides. ft is clear that the scatter of the 
eclipse timings with respect to these ephemerides is far larger 
than can be attributed to measurement error (which can be 
estimated from the scatter in the groups of timings made 
within a few days of each other) . 

The question might be raised as to why we have not 
converted all the timings to Barycentric Julian Dynamical 
Date. The reason is that, as already shown, there are sub- 
stantial residuals, far in excess of measurement error. The 
difference between Barycentric timings and Heliocentric tim- 
ings is a linear shift in the time base (due to the addition of 
leap seconds in the UTC time system), and a periodic error 
of at most 2 s (with roughly the orbital period of Jupiter). 
This periodic error is quite small compared to the residuals 
seen in Fig. 3(a,b) (~ 12 s). We thus prefer to use Heliocen- 
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Figure 3. (a) Fit of a linear ephemeris to the ingress (open cir- 
cles and egress timings (open squares), (b) Fit of a quadratic 
ephemeris to the ingress and egress timings, (c) The mean pe- 
riod estimated from the ingress (open circles) and egress (open 
squares) data. The zeropoint on the vertical axis is at 0.150757525 
d. The points show raw estimates; an interpolated retrospectively 
smoothed estimate passes through the points (e.g. Koen 1996, 
and references therein). For all 3 panels the Julian Date on the 
abscissa is with respect to 2440000.0. 



Table 2. Ephemerides for EC13471-1258 
Parameter 



Epoch 
(JD©) 

Period 
(d) 



Ingress 
2440000+ 

8689.63547 
± 1 



Mid-Eclipse 
2440000+ 

8689.64062 
± 1 



Egress 
2440000+ 

8689.64577 
± 1 



0.150757525 0.150757525 0.150757525 



± 1 



± 1 



± 1 



trie timings for the present and if future observations and 
analysis warrant it, it will be relatively straightforward to 
convert the timings listed in Table 1. 

Although the residuals with respect to the quadratic 
ephemeris are smaller than those of the linear ephemeris 
(as expected from including an additional parameter in the 
model), there is little to choose between the two. Further- 
more, we believe that there is little point in deriving a so- 
phisticated functional fit to the eclipse timings in Table 1: 
there is jitter in the binary period of the system which is 
intrinsic and any functional fit now would likely have no pre- 
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dictive power. Instead, we derive simple linear ephemerides 
for mid-ingress, mid-egress and conjunction as used in: 

JD Ingress = EpQch + ( p oriod x E) 
JD Mid-cclip S c = Epoch + (period x E) 
JD E gross = Epo(;h + ( p eriod x E) 

where E is the cycle number. The parameters are listed in 
Table 2. Errors on the quantities are shown in Table 2: the 
systematic behaviour of the residuals precludes any formal 
estimate of these errors. Instead, these error estimates were 
based on the fact that each eclipse can be timed with an 
accuracy of about 1 s and that the observations span nearly 
24000 cycles. There are, however, systematic departures of 
up to ~ ±12 sec with respect to these ephemerides. 

Koen (1996) proposed a statistical model for sparsely 
observed phases of interest (in the present case eclipses) of 
periodic stars. The observed eclipse timings Tj are assumed 
to be determined by three independent processes: (i) the 
slowly evolving mean period Pj of the binary; (ii) small, ran- 
dom, cycle-to-cycle variations tj in the binary period; and 
(iii) measurement error ej. Assuming that the cycle count 
difference between the j-th and (j — l)-th observation is rij, 

Tj -T j - 1 =n j Pj + ^2 e k + ej=rijPj + Uj + ej (1) 

k=N j _ 1 

where ej and ej are zero mean white noise processes with 
variances and o\, and Nj is the cumulative cycle count 
(i.e. Nj = iVj_i + rij). The specification is completed by 
a model for Pj-. Koen(1996) showed by simulation that the 
random walk assumption 

N 3 

Pj = Pj-i+ % (2) 

k=N j _ 1 

(where r)k is zero mean white noise with variance ajj) works 
very well. It follows that the variances of the three terms on 
the right hand side of (1) are 

vai(njPj) = n|a^ var ({/,-) = Tija\ var(ej) = o\ 

Eqns. (1) and (2) can be solved with the aid of the 
Kalman filtering algorithm; see Koen (1996) for details. The 
optimal solution (according to the Bayes Information Crite- 
rion - see any modern text on time series analysis, or Koen 
1996) has o t = 0, o v = 8.7E - 10 and cr e = 1.05S - 5 (i.e. 
0.9 s), based on the ingress data. The same model is selected 
for the egress data; parameter values are a v = 7AE — 10 and 
o e = 1.36S — 5 (1.2 s). The solutions for Pj are plotted in 
Fig. 3(c). These results show that the period variations are 
not due to cycle to cycle variations but due to a random 
walk in the binary period. Although there is a general trend 
of period shortening in Fig. 3(c) with substantial variations 
about this trend, this behaviour is applicable only to the 
observation period; inferences about the future behaviour of 
the orbital period cannot be drawn from these results. 

2.2 Eclipse analysis 

The white light eclipse light curves (which are in extinction 
corrected counts s" 1 ) were normalized by the count rate 
just prior to ingress and the eclipse ephemeris was used to 
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Figure 4. (a) Mean white light eclipse curve with a model super- 
imposed as described in the text. The ordinate is in units of the 
brightness just prior to eclipse, (b) As for (a) except that only 
the egress section is plotted (on a larger scale). Model fits to the 
egress for i = 75° and i = 90° are superimposed on the data. The 
data and the model for i = 90° have been shifted vertically by 
0.3 for clarity. 

fold all the white light eclipses into 1000 orbital phase bins 
between phase -0.05 and 0.05. Small corrections were used 
to the eclipse timings to reduce the intrinsic jitter and pre- 
cisely "align" the eclipses. The result is shown in Fig. 4(a). 
The resolution of 0.0001 in phase corresponds to ~ 1.3 s, 
comparable to the original time resolution of each individ- 
ual eclipse. The reduction in scatter compared to the eclipse 
plotted in Fig. 1 is obvious. Analysis of the light curves 
of totally eclipsing binaries can, in principle, provide iii /a, 
Rila and i where, in the present context, Ri is the radius 
of the white dwarf, R2 is the radius of the M dwarf, a is the 
orbital separation and i is the inclination (throughout this 
paper, subscripts 1 and 2 shall refer to the white dwarf and 
M dwarf respectively). However, Ritter & Schroder (1979: 
RS79) have shown that because eclipses of the kind seen in 
Fig. 4(a) have only two well-defined properties, the times 
from first to fourth and second to third contacts, the anal- 
ysis of such eclipses cannot provide a unique solution to all 
three quantities. For the sake of brevity, we do not repeat 
RS79's arguments here and refer the reader to their paper. 
We do, however, follow their analysis closely and we fit the 
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mean eclipse in Fig. 4(a) (solid line) with a slightly different 
analytic function to that which they derived: 

I it) = A + A 1 (t-to)+d.f{Ri/a,R 2 /a,i,t-to) 

f = 1, for x < —p 

f = 0.5 + l/7rp 2 [p 2 sin _1 (s/p) 

. , 2 2-.1/2 ,/ 2 2\/-, 2 . 2x11/2 

+s(p -s ) -[(p -s)(l-p +s)]' 
+sin _1 (p 2 — s 2 ) 1 ^ 2 ] for —p<x<p 

f = 0, for a; > p 

s = (2x - x 2 - p 2 )/(2(l - x)) 
p = R1/R2 

x = 1- (4ir 2 /P 2 (t-t ) 2 + cos 2 i) 1/2 /(R 2 / a) 

In the above equations, I(t) describes the brightness of the 
system near and during eclipse. Ao is the brightness dur- 
ing eclipse and d is the depth of the eclipse. A\ allows for 
a linear variation in brightness during the observation. / 
is a function that varies between and 1 according to the 
value of x: x < — p outside eclipse and x > p during total- 
ity. P is the orbital period and to is the time of mid eclipse. 
The indeterminacy of these equations was resolved by choos- 
ing a value for the inclination angle i. The function I was 
then fitted to the light curve using a nonlinear least squares 
technique, optimizing 6 parameters: Ao, Ai, d, to, Ri/a and 
R 2 /a. These calculations were performed for a range of i 
from 70° to 90° in steps of 1°. The parameters Ao,Ai,to 
and d did not vary with i. Their values were found to be: 
A = 0.247 ± 0.001, d = 0.755 ± O.OOli, ; Ai = 0.090 ± 0.001 
(in units of the out-of-eclipse brightness) and to = 0.0 s. 
The positive slope, Ai, during totality is significant and will 
be discussed further below. The most interesting quantities 
from the eclipse analysis are the stellar radii. These are listed 
in Table 3 as a function of i. The model fit for inclination 
75° (an arbitrary choice) is shown as the solid line in Fig. 
4(a). 

Fig. 4(b) shows the mean eclipse egress on a larger scale 
with the models for i = 75° and i = 90° superimposed. 
Close inspection shows that the models are in satisfactory 
agreement with the data and, as pointed out by RS79, es- 
sentially identical. Even with the high signal-to-noise in the 
mean eclipse curve, the data are incapable of distinguish- 
ing between the two models. The models for 75° and 90° 
were divided by each other to show the subtle difference be- 
tween the models. The data were also divided by each model 
in turn to see if any such trend could be discerned. Again, 
there was no basis from these operations on which to fix the 
inclination. 

The traditional eclipse parameters: the time from first 
to fourth contact (0m), the time from second to third con- 
tact (023), and the time from mid-ingress to mid-egress 
of the primary (^1/2) (e.g. see equations 2-4 of Robinson, 
Nather & Patterson 1978), were determined to be: 

0i4 = 0.07246 ± 0.00010 (944 ± 1.3s) 
023 = 0.064 08 ± 0.00005 (835 ± 1.0s) 
0i/2 = 0.06832 ± 0.00001 (889.9 ± 0.2s) 



Table 3. Eclipse parameters 



i 


R 2 /a 


Ri/a 


70 


0.404 





.00695 


71 


0.390 





.00720 


72 


0.376 





.00746 


73 


0.363 





.00775 


74 


0.349 





.00805 


75 


0.336 





.00837 


76 


0.323 





.00870 


77 


0.311 





.00906 


78 


0.299 





.00944 


79 


0.287 





.00983 


80 


0.276 





.01023 


81 


0.266 





.01064 


82 


0.256 





.01106 


83 


0.247 





.01148 


84 


0.239 





.01187 


85 


0.232 





.01225 


86 


0.226 





.01258 


87 


0.221 





.01286 


88 


0.217 





.01307 


89 


0.215 





.01320 


90 


0.215 





.01325 



All good quality individual runs were also fitted with 
the model. No significant variations in eclipse widths were 
found. There was some variation in the eclipse depths but 
given the uncontrolled photometry (different telescopes and 
detectors as well as no observations of standard stars) , these 
modest variations cannot be attributed to the star. However, 
the significant positive slope seen in the model of the mean 
eclipse was found in many of the individual runs (as expected 
as the mean eclipse was derived from the individual runs). 
This positive slope is certainly not of instrumental origin but 
arises from the binary. In roughly fifty per cent of high speed 
photometry runs, a positive slope (of variable size) during 
eclipse was detected; no negative slope was ever detected; in 
the remaining cases, the slope during eclipse was consistent 
with zero. We discuss this point further in the section on 
multicolour photometry. 

2.3 Flares 

The ellipsoidal modulation shown by the companion M 
dwarf (Fig.l) indicates that it is tidally locked to the bi- 
nary period and that therefore it is certainly rotating very 
rapidly. However, for M dwarfs (and later types), it is not at 
all clear that chromospheric activity is correlated with rota- 
tion (see, for example, the recent review by Hawley, Reid & 
Gizis 2000). The presence of Ha emission - usually regarded 
as the primary indicator of magnetic activity - is not affected 
by heating from the very close white dwarf (see later). Nev- 
ertheless, if we assume that the M dwarf is similar to GJ 
486 (as indicated by the colour comparison given in section 
3.2 and later discussion), then the implied spectral type of 
M3.5 means that the M dwarf is completely convective and 
in a regime where at least 20 per cent of all stars show Ha 
emission and an even greater fraction exhibit flaring activ- 
ity over a large range of energies (see Hawley, Reid & Gizis 
2000; especially their Figs. 2 and 3). The situation with the 
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Table 4. Photoelectric UBV photometry 

JD© Vtot B — Vtot U — B tot B tot t/tot 

2440000+ 

1992 Mar 04/05 



8686.634 


14.42 


0.63 


-0.54 


15 


.05 


14.51 






1992 Mar 06/07 








8688.544 


14.37 


0.63 


-0.54 


15 


.00 


14.46 


8688.570 


14.47 


0.58 


-0.58 


15 


.05 


14.47 


8688.596 


14.43 


0.56 


-0.54 


15 


.00 


14.45 






1992 Mar 08/09 








8690.545 


15.21 


1.54 


0.36 


16 


.75 


17.10 


8690.564 


14.41 


0.59 


-0.54 


15 


.01 


14.47 






1992 Mar 09/10 








8691.458 


7.15* 


0.42 


0.07 








8691.473 


14.52 


0.59 


-0.53 


15 


.11 


14.58 


8691.481 


14.36 


0.65 


-0.52 


15 


.00 


14.49 


8691.485 


14.36 


0.65 


-0.52 


15 


.01 


14.49 


8691.493 


7.13* 


0.42 


0.07 








8691.497 


14.34 


0.65 


-0.54 


14 


.99 


14.45 


8691.502 


14.39 


0.60 


-0.51 


14 


.98 


14.47 


8691.511 


14.45 


0.56 


-0.55 


15 


.01 


14.47 


8691.516 


14.46 


0.56 


-0.66 


15 


.09 


14.36 


8691.524 


7.13* 


0.42 


0.07 








8691.528 


14.49 


0.55 


-0.60 


15 


.04 


14.44 


8691.534 


14.48 


0.55 


-0.57 


15 


.03 


14.45 


8691.539 


14.44 


0.57 


-0.57 


15 


.01 


14.44 






1992 May 29/30 








8772.222 


7.15* 


0.40 


0.06 








8772.229 


14.45 


0.58 


-0.52 


15 


.03 


14.50 


8772.235 


14.42 


0.60 


-0.52 


15 


.02 


14.51 


8772.241 


14.47 


0.58 


-0.49 


15 


.05 


14.55 


8772.271 


7.14* 


0.41 


0.05 








8772.282 


14.42 


0.59 


-0.50 


15 


.01 


14.51 


8772.289 


14.38 


0.62 


-0.53 


15 


.00 


14.47 


8772.295 


14.39 


0.60 


-0.53 


14 


.98 


14.45 


8772.301 


7.15* 


0.41 


0.06 








8772.304 


14.43 


0.57 


-0.54 


15 


.00 


14.47 


8772.310 


14.44 


0.59 


-0.56 


15 


.03 


14.47 


8772.316 


14.45 


0.59 


-0.56 


15 


.03 


14.47 


8772.321 


7.15* 


0.41 


0.06 








8772.325 


14.48 


0.47 


-0.68 


14 


.95 


14.27 


8772.331 


14.48 


0.53 


-0.63 


15 


.02 


14.38 


8772.337 


14.50 


0.52 


-0.60 


15 


.03 


14.43 


8772.343 


7.15* 


0.41 


0.06 








8772.346 


14.46 


0.58 


-0.60 


15 


.04 


14.44 


8772.352 


14.46 


0.56 


-0.54 


15 


.02 


14.48 


8772.358 


14.44 


0.57 


-0.56 


15 


.01 


14.45 


8772.364 


7.15* 


0.41 


0.07 








8772.367 


14.39 


0.62 


-0.51 


15 


.01 


14.51 


8772.374 


14.46 


0.57 


-0.49 


15 


.03 


14.54 


8772.380 


14.45 


0.61 


-0.51 


15 


.06 


14.55 


8772.385 


7.18* 


0.42 


0.06 








8772.388 


14.51 


0.58 


-0.57 


15 


.08 


14.52 


8772.394 


14.54 


0.53 


-0.51 


15 


.07 


14.56 


8772.421 


14.50 


0.58 


-0.53 


15 


.08 


14.55 


8772.426 


7.17* 


0.42 


0.08 








8772.429 


14.45 


0.58 


-0.52 


15 


.03 


14.51 


8772.435 


14.40 


0.63 


-0.51 


15 


.04 


14.53 


8772.440 


14.40 


0.61 


-0.54 


15 


.01 


14.47 


8772.445 


7.17* 


0.42 


0.07 








8772.449 


14.39 


0.61 


-0.53 


15 


.00 


14.46 


8772.454 


14.45 


0.57 


-0.51 


15 


.02 


14.51 


8772.460 


14.42 


0.64 


-0.58 


15 


.06 


14.48 


8772.466 


7.19* 


0.42 


0.08 









Table 4 - continued 



JD 


Vtot 




B - Vtot 


U - Btot 


Btot 


ftot 


2440000+ 




















1997 Apr 29/30 






10568.358 


14.39 




0.63 


-0.55 


14.99 


14.44 


10568.366 


14.40 




0.63 


-0.51 


15.03 


14.51 


10568.371 


14.43 




0.60 


-0.53 


15.03 


14.50 


10568.379 


15.24 




1.55 


0.88 


16.79 


17.67 


10568.384 


15.22 




1.52 


0.70 


16.74 


17.43 


10568.389 


14.43 




0.59 


-0.53 


15.02 


14.49 


10568.394 


14.40 




0.64 


-0.53 


15.04 


14.52 








1997 May 03/04 






10572.285 


14.39 






n ^9 
-u.oz 


15.04 


14.52 




14.39 




0.63 


-0.53 


lO.UO 


14 49 


10572.298 


15.22 




1.56 


0.98 


16.78 


17.76 


10572.302 


15.21 




i fii 


n on 


16.82 


17.72 


iuo i z.ouy 


14.42 




0.63 


-0.51 


LO.UO 


14 54 


10572.313 


14.39 




0.65 


-0.50 


15.04 


14.54 


10572.320 


14.38 




0.65 


-0.49 


15.02 


14.54 


10572.324 


14.39 




0.66 


-0.52 


15.04 


14.52 




14.43 




0.59 


-0.54 


15.03 


14.47 




B — Vvvd 


U — Bwd 


%D 


UwD 




ouyu. o<±o 


15.15 




0.09 


-0.68 


15 24 


1 A ^fi 


10568.379 


15.13 




0.13 


-0.70 


15.26 


14.56 


1UOOO. OOQ 


15.15 




0.12 


-0.70 


1 ^ 97 


1 A ^7 


10572.298 


15.11 




0.17 


-0.71 


15.28 


14.57 


10572.302 


15.12 




0.15 


-0.70 


15.27 


14.57 




15.13 




0.13 


-0.70 


15.26 


14.57 




Vsec 




B - V Scc 


U - B Sec 


-Bscc 


^Soc 


8690.545 


15.21 




1.54 


0.36 : 


16.75 


17.10 : 


10568.379 


15.24 




1.55 


0.88 : 


16.79 


17.67 : 


10568.384 


15.22 




1.52 


0.70 : 


16.74 


17.43 : 


10572.298 


15.22 




1.56 


0.98 : 


16.78 


17.76 : 


10572.302 


15.21 




1.61 


0.90 : 


16.82 


17.72 : 




15.22 




1.56 


0.76 : 


16.78 


17.54 : 


Table 5. Single filter time series photometry 








Run 


Date Start of Run 


Length 


Filter 




Name 




JD© 


(hi) 










2440000+ 








S5502 


29 


5 92 8772 


.245553 


0.5 


V 




S5503 


29 


5 92 8772 


,398321 


0.5 


U 




S5601 


28 


2 93 9047 


.481823 


3.8 


V 




S5603 


01 


3 93 9048 


.433285 


0.5 


B 




S5603a 


01 


3 93 9048 


.454120 


4.4 


V 




ck22 


14 


4 93 9092 


.426764 


3.8 


V 




ck25 


17 


4 93 9095 


.414611 


2.4 


V 




ck26 


18 


4 93 9096 


.339396 


5.4 


V 




m9460d 


17 


4 94 9460 


.406300 


4.1 


V 




S5758 


11 


6 94 9515 


.325139 


3.2 


V 
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ck70 (Top) S5853 (Bottom) 

(Sec) 

500 1000 1500 



1500 
1000 
500 





Figure 5. (Upper) Run ck70 showing a substantial and double- 
peaked flare. The ordinate is in extinction corrected counts s _1 . 
(Lower) Run S5853 also shows a single-peaked flare of lower am- 
plitude. 



cool star in EC 13471-1258 may be complicated by the fact 
that it will have certainly passed through a "common en- 
velope" configuration and may have undergone substantial 
mass loss through Roche lobe overflow (see the later discus- 
sion). Nonetheless, it might still be expected to find flares 
of the kind seen in chromospherically active M dwarfs. 

Unequivocal evidence for flares is seen in the two light 
curves (ck70 and S5853) plotted in Fig. 5. Both flares began 
in eclipse, thus identifying the red star as their origin. The 
flare in ck70 is ~ 3 mag brighter at maximum and the profile 
is double peaked (the third peak is the egress of the white 
dwarf). The flare in S5853 is only 1.5 mag in size. The ck70 
light curve continued beyond what is plotted in Fig. 5 with 
the brightness of the system gradually returning to the pre- 
eclipse brightness level. The rise time of the lower amplitude 
flare is a few seconds whereas in ck70 it was about a minute. 

We have observed two flares out of a total of ~ 40 hours 
of high speed photometry. Two additional flares were ob- 
served in multicolour photometry (see next section) cover- 
ing 21 hours. Flaring therefore occurs with a mean rate of 
0.067 per hour. We have also seen evidence for flares in spec- 
troscopy (as did Kawka et al. 2002), although it is hard to 
quantify the rate. 



3 MULTICOLOUR PHOTOMETRY 

We now consider multicolour photometric data in order to 
learn about variations due to the M dwarf and the rela- 
tive contributions of each component star at different wave- 
lengths. 



> 



m 



14.3 =- 

14.4 =- 

14.5 1- 

14.6 f- 

14.9 j- 

15 =- 

15.1 =- 

15.2 =- 

14.3 I- 

14.4 =- 

14.5 1- 

14.6 =- 



8772.2 8772.3 8772.4 8772.5 
HJD 

Figure 6. UBV photometry on 29/30 May 1992 (JD 2448772). 
Note the flare in the U and B light curves. 



3.1 Photoelectric photometry 

In order to assist any future modelling of the binary, we list 
in Table 4 all the photoelectric UBV measurements. They 
were made with blue sensitive photomultiplier tubes. The 
data were reduced by correcting for atmospheric extinction 
and transforming to the standard system using observations 
of several tens of E region standard stars (Menzies et al. 
1989). Standard stars were also observed regularly during 
the nights listed in Table 4. In addition, on the nights of 1992 
Mar 9/10 and May 29/30, a ~ 7th mag local comparison star 
HD120543, about 10 arcmin distant from EC13471-1258, 
was also observed. Its magnitude and colours are shown by 
the data in Table 4 with a superscript asterisk (*) along- 
side. These data are included so the reader can judge the 
quality of the photometry. Note that towards the end of the 
night 1992 May 29/30 the comparison star faded by a few 
per cent. The EC13471-1258 data were corrected by these 
differences, on the assumption that the comparison star is 
constant. Even if this is not so, the corrections were at most 
0.03 mag in V and no correction at all in the colours. We 
estimate the uncertainty on the photometry is ±0.02 mag. 

Almost all the data in Table 4 are labelled with the 
subscript "tot" to indicate the sum of the light from both 
component stars. Five data points were obtained during to- 
tal eclipse of the white dwarf, allowing a separation of the 
contributions of each star. These separate contributions are 
shown at the bottom of Table 4 with subscripted labels 
"WD" and "Sec". Bold numbers in Table 4 are means of 
the respective quantities. In the case of the "tot" data, the 
eclipse points were excluded. The eclipse depths in V, B and 
U are 56, 81 and 93 per cent respectively. Only 7 per cent of 
the U light originates from the M dwarf, so its U magnitude 
and U — B colour are not well determined; this is denoted 
in Table 4 by colons. 

Fig. 6 shows the data obtained on 1992 May 29/30 (JD 
2448772) excluding the eclipses which were observed in the 
single filter runs S5502 and S5503 (Table 5). The V light 
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14.4 
14.5 
14.6 

14.9 
15 
15.1 
15.2 

14.3 
14.4 
14.5 
14.6 



X- V." , ** > ; 



.2 .4 .6 .8 1 
Orbital Phase 

Figure 7. All UBV photometry listed in Table 4 folded on the 
orbital period. Eclipse points have been excluded. 



curve shows the ellipsoidal variation (double humped struc- 
ture) already pointed out in Fig. 1 (upper). In this case, 
however, the variation is more erratic (that this is real can 
be judged from the photometry of the comparison star: Ta- 
ble 4). The ellipsoidal variation is less obvious in the B and 
U light curves. Instead, the U light curve shows a flare which 
is less evident in B and undetected at V. Such behaviour is 
typical of flare stars where increases in optical light are most 
obvious due to the Balmer continuum and lines going into 
emission. Fig. 7 shows all the UBV data in Table 4 folded 
on the orbital period. As in Fig. 6, the ellipsoidal variation 
is obvious only in the V light curve. The U band light curve 
seems to reach a maximum at phase 0.5 (even after allowing 
for the occurrence of a flare close to this phase). The appar- 
ent increase in scatter in all the magnitudes close to eclipse 
is real: the effect is subtle but can be best appreciated by 
comparing the light curves at orbital phase 0.5 with phases 
close to 0. 

Table 5 lists single filter observations made on the 
SAAO 0.75-m and 1-m telescopes with photomultiplier de- 
tectors. With the exception of the first two runs, S5502 and 
S5503, which used 1 s integrations, all the other data were 
obtained with 10 s integrations. The photomultiplier data 
were reduced in the same manner as the white light data 
except that coefficients of 0.30 and 0.54 mag/airmass were 
used for atmospheric extinction correction of the B and U 
data respectively (these are close to the standard extinction 
coefficients for the site). 

Fig. 8 (upper) shows the V data from Table 5. A number 
of features are apparent: 

• Taking into account that 56 per cent of the light in 
the V band data is due to the white dwarf, the ellipsoidal 
variation is of large amplitude: a sinusoid with period equal 
to half the orbital period was fitted to the data in Fig. 8. The 
mean amplitude (peak to peak) of the ellipsoidal variation 
was found to be 22 per cent of the mean brightness of the 
M dwarf. 



600 
400 
200 




EC13471-1258 V Data 

PHASE 

3.5 1 1.5 

—I 1 1 1 1 1 1 1 1 r 



M ,„„,. M ********* 

— I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 V 

to* 



— I 1 1 H 



H I I 1 I I 1 I I 1 I I 1 1 I 1 h 



600 



JD 9092+ 
0.45 0.455 



(Days) 
0.46 0.465 



Figure 8. (Upper) V filter time series photometry folded on the 
orbital period. The runs are listed in Table 5 and are from the 
top: S5601, S5603a, ck22, ck25, ck26, m9460d and S5758. (Lower) 
Run ck22 close to eclipse and plotted on a larger scale. Note the 
two pre-eclipse dips. 



• The minimum of the ellipsoidal variation is not centred 
on eclipse. This is most noticeable in runs ck22 and ck25 (3rd 
and 4th from the top). This gives rise to the positive slope 
during total eclipse already discussed in the white light pho- 
tometry. The minimum closest to eclipse in the least squares 
fit mentioned above occurred at orbital phase —0.02340.001. 

• Although of low amplitude, there appears to be rapid, 
erratic variations superimposed on the overall ellipsoidal 
modulation. This is most apparent again in ck22 and ck25. 
For example, in ck22 (third panel from top) at orbital phase 
~ 1.2 there are rapid short time scale variations which we 
believe are not due to atmospheric transparency variations 
or of instrumental origin. 

• ck22 shows two pre-eclipse dips, plotted in Fig. 8 
(lower). We are convinced that these are real and not of 
instrumental origin: the dip minima are at the level of the 
eclipse and not the sky background so the effect cannot be 
explained by the star going out of the photometer's aperture. 
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Table 6. CCD VRI photometry of EC13471-1258 

JD Q V JD Q R JD q I 

2440000+ 2440000+ 2440000+ 



9161 


.351 


14. 


38 


9161 


.352 


13 


.67 


9161 


.353 


12 


.24 


9161 


.354 


14. 


38 


9161 


.355 


13 


.66 


9161 


.359 


12 


.32 


9161 


.357 


15. 


20 


9161 


.358 


13 


.99 


9161 


.363 


12 


.31 


9161 


.361 


15. 


22 


9161 


.362 


13 


.98 


9161 


.366 


12 


.28 


9161 


.364 


15. 


17 


9161 


.365 


13 


.97 


9161 


.378 


12 


.16 


9161 


.367 


14. 


36 


9161 


.377 


13 


.58 


9161 


.381 


12 


.15 


9161 


.375 


14. 


32 


9161 


.380 


13 


.56 


9161 


.384 


12 


.12 


9161 


.379 


14. 


32 


9161 


.383 


13 


.54 


9161 


.388 


12 


.09 


9161 


.382 


14. 


29 


9161 


.387 


13 


.51 


9161 


.391 


12 


.08 


9161 


.385 


14. 


28 


9161 


.390 


13 


.49 


9161 


.394 


12 


.06 


9161 


.389 


14. 


28 


9161 


.393 


13 


.47 


9161 


.398 


12 


.05 


9161 


.392 


14. 


25 


9161 


.397 


13 


.46 


9161 


.401 


12 


.04 


9161 


.396 


14. 


25 


9161 


.400 


13 


.46 


9161 


.404 


12 


.06 


9161 


.399 


14. 


25 


9161 


.403 


13 


.47 


9161 


.408 


12 


.09 


9161 


.402 


14. 


24 


9161 


.407 


13 


.47 


9161 


.411 


12 


.11 


9161 


.406 


14. 


26 


9161 


.410 


13 


.51 


9161 


.415 


12 


.13 


9161 


.409 


14. 


.27 


9161 


.414 


13 


.53 


9162 


.339 


12 


.33 


9161 


.412 


14. 


25 


9162 


.338 


13 


.73 


9162 


.343 


12 


.33 


9162 


.337 


14. 


41 


9162 


.342 


13 


.73 


9162 


.346 


12 


.33 


9162 


.340 


14. 


42 


9162 


,345 


13 


.73 


9162 


.349 


12 


.31 


9162 


.344 


14. 


40 


9162 


.348 


13 


.72 


9162 


.353 


12 


.27 


9162 


.347 


14. 


39 


9162 


.352 


13 


.69 


9162 


.356 


12 


.24 


9162 


.350 


14. 


38 


9162 


.355 


13 


.67 


9162 


.361 


12 


.19 


9162 
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Table 7. VRI magnitudes and colours of the component stars 




^ojP 14.28 



JD 2449163 




JD 2449164 «*, I 3 - 10 
^ ^ 13.53 

r m rffD 14.26 
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Xxxxx*x»xx*<xx xj< x x xxxxxxxxxx^xx xX 
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.5 1 1.5 2 

Orbital Phase 

Figure 9. CCD VRI photometry on HJD 2449161, 2449162, 
2449163 and 2449164. Asterisks indicate I magnitudes, filled cir- 
cles R magnitudes and open circles V magnitudes. Orbital phase 
with respect to the ephemeris in Table 2 is indicated on the ab- 
scissa. Ordinate carets are spaced at 0.5 mag intervals. The zero 
point for each filter subset is different but can be deduced from 
the magnitude indicated next to the last point in each filter sub- 
set. The lowest curve (crosses) shows the V — I colour on HJD 
2449163. Note the decaying flare at the start of the data on HJD 
2449163. 



Other tests were done by two of us (CK and DOD) during 
this observing run to search for any possible atmospheric 
or instrumental origin: no evidence was found. We accept, 
however, that there was no comparison star observations to 
prove the reality of these dips and they were observed only 
in ck22. 



3.2 CCD photometry 

CCD VRI photometry was obtained over the four nights 
1993 June 22-25 (JD 2449161-2449164) on the 1-m telescope 
using the UCL camera with the RCA chip. Repeated cycles 
of V, R and / exposures, with respective integration times 
of 100, 80 and 60 s were obtained. Prefiash and flat field 
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calibration frames were also obtained (the latter from the 
twilight sky in photometric conditions). 

Observations on the first three nights were hampered 
by cirrus so differential magnitudes for EC13471-1258 with 
respect to two local comparison stars on the frame were ob- 
tained. The fourth night was clear: E-region standard stars 
(Menzies et al. 1989) were observed and enabled transforma- 
tion of the instrumental magnitudes of the comparison stars 
and EC13471-1258 to the V(RI)c system. We estimate the 
uncertainty on the photometry is ±0.03 mag. About 10.6 hr 
of useful data were obtained. Once again, in order to assist 
modelling of the binary, the VRI photometry are listed in 
full in Table 6. 

The results are plotted against orbital phase (from Ta- 
ble 2) in Fig. 9. Most obvious is the previously mentioned 
ellipsoidal variation and the eclipses. The eclipse depths de- 
crease with wavelength: in VRI the depths are respectively 
53, 26 and 7 per cent. Note the decay of a flare at the start 
of the run on JD 2449163: the decline lasted ~ 1 hr, similar 
to the flare in Fig. 5 (upper). 

The ellipsoidal variation was determined by fitting a 
sinusoid with half the orbital period to the data with the 
eclipses and the flare excised. After correcting for the white 
dwarf contribution, the amplitude (peak to peak) of the fit- 
ted sinusoid in V, R and I was respectively 27, 25 and 22 
per cent. The V amplitude from the photoelectric photom- 
etry was 22 per cent. However, there are short term varia- 
tions in the ellipsoidal variation: for example, in Fig. 9, on 
JD2449163, the first hump after the flare is of lower ampli- 
tude than the second hump. 

Four eclipses were observed during the VRI CCD pho- 
tometry and this allowed separation of the magnitudes and 
colours of the component stars; the results are shown in Ta- 
ble 7 along with the means of the four determinations. The / 
and V — I results for the white dwarf are uncertain because 
it contributes less than 10 per cent in the / pass band. Note 
that the V magnitude of the white dwarf as determined from 
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Figure 10. Upper curve: sum of G140L spectra of EC13471-1258 
on 1999 Aug 28. The individual spectra were velocity shifted to 
the white dwarf rest frame prior to summation. Various absorp- 
tion features are marked. Middle curve: best fit model with T e jf 
14220 K and log g of 8.34. The model has been offset downwards 
by 0.8 X 10~ 14 with respect to the observations for clarity. Lower 
curve: residuals of the fit in the sense of (data - model). The 
residuals have also been shifted downwards 1.0 X 10~ 14 for clar- 
ity. Further details are given in the text. The abscissa is in A and 
the ordinate is in erg s _1 cm~ 2 A -1 . 

the CCD data is 0.1 mag brighter than determined by the 
UBV data (Table 4). We felt unwilling to force agreement 
between the two: the difference can be regarded as an es- 
timate of the total external error of the measurements. We 
prefer the UBV determination as the process of multicolour 
photometry was much better calibrated for the photoelectric 
data compared to the CCD data. 

Also shown in Table 7 are our measurements, using the 
same equipment, of the colours of three nearby M dwarfs: 
GJ 486, 514 and 551. A full discussion of the secondary star 
will be given later but the results in Table 7 show that the 
M dwarf in EC13471-1258 is similar, at least as far as broad 
band colours are concerned, to GJ 486. 



4 SPECTROSCOPY 

4.1 Ultraviolet Spectroscopy 

In order to determine the effective temperature and grav- 
ity of the white dwarf, as well as measure its radial velocity 
curve, HST/STIS observations of EC13471-1258 were ob- 
tained during a single visit lasting 3 orbits on 28 August 
1999. All observations were made through the 52" x 0.2" 
aperture. In the first orbit, three spectra of 300 sec exposure 
time were obtained with the G140L grating with resolution 
of ~ 1.5 A and wavelength coverage from 1150 to 1700 A. 
The remainder of the first orbit was used to obtain two 300- 
sec spectra with the G230L grating covering the wavelength 
range 1600-3200 A with 3.3 A resolution. In each of the sec- 
ond and third orbits, six 300-sec spectra were obtained with 
the G140L grating in the same manner as in the first orbit. 

With the exception of CI at 1660 A, and the Mgll dou- 
blet at 2800 A, the G230L data arc featureless. In contrast, 
the uppermost curve in Fig. 10 shows the sum of the fifteen 
G140L spectra. Prior to summation, the spectra were shifted 



into the rest frame of the white dwarf using the radial veloc- 
ity solution determined below. This procedure avoided the 
small but detectable smearing of the spectrum due to the 
white dwarf's orbital motion. The summed spectrum shows 
a continuum rising to short wavelengths but interrupted by 
very strong Lya absorption (with narrow geocoronal emis- 
sion in the core). This, along with the quasimolecular 
feature at 1400 A (Allard & Koester 1992) indicates that 
the white dwarf is cool (well under 20000 K). In addition, 
there are narrow metal absorption lines, mostly due to CI 
and Sill. The middle curve in Fig. 10 shows our best-fitting 
model as described in the next subsection; it has been dis- 
placed downwards by 0.8 x 10 -14 for clarity. Although a good 
match to the observations is evident, the fit is not perfect. 
We defer detailed discussion of the inadequacies of the fit 
until the end of the next subsection. 

4-1.1 Modelling The Spectrum of the White Dwarf 

The effective temperature, surface gravity, chemical abun- 
dance and rotational velocity of the white dwarf were de- 
termined by the classical technique of model atmosphere 
calculation and fitting of the synthetic spectrum to the ob- 
servational data. 

The theoretical models we employed in these calcula- 
tions use the same procedures and programs developed over 
many years and applied to all spectral classes of white dwarfs 
by the Kiel group. The codes and the input physics are very 
similar to the description in Finley et al. (1997). A minor 
change compared to most applications is the chemical com- 
position. As the UV spectrum shows metal lines, we have 
computed a new grid of hydrogen-rich models with metals 
in relative solar abundance. 

In order to limit the size of the grid of models, some 
initial estimates of the atmospheric parameters were made 
by trial and error. As mentioned before, the presence of the 
quasi-molecular feature at 1400 A, as well as the absence 
or marginal detection of the Hb feature at 1600 A enabled us 
to set upper and lower bounds to the required temperature 
range. An initial estimate of the line broadening was deter- 
mined by convolving the theoretical spectra with Gaussians 
of varying FWHM. A value of ~ 3 A gave a good initial 
fit. With this in hand, it rapidly became apparent that a 
solar metal abundance resulted in metal lines far too strong 
compared to the observations; on the other hand, an abun- 
dance of 1/100 solar produced lines of insufficient strength. 
The grid that was computed, therefore, had effective tem- 
peratures in the range 13500 to 16000 K in steps of 250 K, 
surface gravities, log g, in the range 7.5 to 8.5 in steps of 
0.25, and metal abundances of 1/3, 1/10, 1/30, and 1/100 
solar. 

The best model fit within this grid was determined with 
a Levenberg-Marquardt \ 2 algorithm (Press et al. 1992). 
Our version of the method is further described in Homeier 
et al. (1998). In order to restrict the range of possible mod- 
els even further, we constrained the fit to be consistent with 
the V magnitude of the white dwarf, derived from both the 
photoelectric and CCD photometry during eclipse (Tables 
4 and 7). We used the mean of the two measurement tech- 
niques V = 15.08 and assigned an error estimate of ±0.05 
mag. The magnitude was converted to flux using the zero 
point in Bessell (1979). 
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The atmospheric parameters so derived were T c ff = 
14220 K and log g = 8.34, mainly determined by the slope 
of the UV flux, the optical V magnitude, and the broad 
features of the Lyman a line and its 1400 A satellite. The 
weak metal lines do not significantly influence this determi- 
nation; we used the grid with 1/30 solar abundances, except 
for Si, which was decreased by an additional factor of ap- 
proximately 2. We discuss uncertainties in these estimates 
in the next subsection. 

The modelling process also showed that the widths of 
the metal absorption lines could not be explained by the 
~ 1.5 A instrumental resolution alone. It is, of course, possi- 
ble that the instrumental resolution is somewhat worse than 
nominal. However, if present, any degradation can only be 
modest as the radial velocity accuracy of the results in the 
next subsection would not have been possible if the instru- 
mental resolution was severely degraded. The most obvious 
explanation for the additional broadening required is rota- 
tional motion of the white dwarf. This was investigated by 
first broadening the best model with Doppler profiles of 0, 
300, 400, 500, 600, 700 and 800 km s _1 . 

Rotational broadening of spectral lines is described by 
a convolution of the intrinsic profile with the broadening 
function 



A(x) 



l + |/3 



(see e.g. Unsold 1968). Here, the variable x is the distance 
from the line center in units of the rotational broadening 



AA 



with AA r0 i = A 



V sini 



AAroi C 

P is the limb darkening coefficient with the angle-dependent 
intensity written as 

/(/0= 7(0) (1 + /?/*) 

with fi = cos 8, and 9 the angle between the line-of-sight and 
the normal to the surface. 

Our code is able to calculate the intensity for different 
angles from the final model, in addition to the emerging flux. 
Such a calculation could be used to estimate (3. There are 
two difficulties with this direct approach: first, the program 
cannot calculate the intensity directly at the limb fj, = 0, 
and second, the limb darkening turns out to be strongly 
non-linear in the UV part of the spectrum. We have there- 
fore estimated (3 by using the outermost point calculated on 
the disk (~ 0.04) as an approximation for the limb value and 
calculating the slope from the intensity at this point and at 
the center of the disk (/i = 1). The values for P decrease 
from ~ 33 near 1150A, to 7 at 1800A and 0.4 in the red 
at 8000 A, with stronger variations through the profiles of 
stronger lines. The classical value for Eddington limb dark- 
ening for a source function increasing linearly with optical 
depth is 1.5. Fortunately, the rotational profile (for an in- 
finitely sharp line) does not change very much for /3 changing 
from to 00, as can be seen from Fig. 168 in Unsold (1968). 
This was confirmed by our own tests, comparing a convolu- 
tion with p = 1 with one with the value appropriate for the 
observed STIS spectral range of P = 15. 

We then convolved each of these spectra with a Gaus- 
sian of FWHM equal to the instrumental resolution of ~1.5 
A. The resulting spectra were both subtracted and divided 



into the observational data and the resulting difference and 
quotient spectra examined (a more formal procedure is not 
justified given the difficulties in detailed matching discussed 
in the next subsection) . If the model lines were too narrow, 
these spectra contained residual absorption; if too deep they 
showed features resembling a sine function (sin x/x). Our 
best estimate of Vi sin i is 400 ± 100 km s~\ 

4-1.2 Error Estimates 

Of equal importance to the atmospheric parameters of the 
best-fitting model is an estimation of the uncertainties, es- 
pecially for the effective temperature and surface gravity. 
Accordingly, this issue was investigated in some detail. 

The residuals between the best-fitting model derived 
in the previous subsection and the summed G140L data 
(shifted to remove the orbital motion of the white dwarf) 
are shown as the lowest curve in Fig. 10. The residuals have 
been shifted down by 1.0 x 10~ 14 for clarity, and are char- 
acterized by: 

• Low frequency trends of up to several per cent on scales 
of tens or hundreds of A, most obvious at 1400 A and at the 
red end of the spectral range shown. 

• High frequency features associated with the metal ab- 
sorption lines. 

We attribute the low frequency trends to a combination 
of (i) an imperfect fit of the quasimolecular Hj feature at 
1400 A. Despite the best available input physics, we were 
unable to achieve a better fit than shown in Fig. 10; and 
(ii) imperfections in the STIS flux calibration, especially at 
the red end of the G140L spectral range. (Note that we do 
not believe that the rise in flux at the red end is due to the 
quasimolecular H2 feature at 1600 A). 

In regard to the poor fits to the metal absorption lines, 
the worst case involves the CI and II lines at A 1330 and 
1335 A where our model has about equal strength for these 
features. In the data, however, the CI A1330 A line is much 
weaker. On the other hand CI A1657 A is stronger in the 
data than in our model. Despite checks of gf values for these 
features, we could find no reason to change the theoretical gf 
values so these disagreements remain unexplained. The Sill 
lines at 1530 A were also found to be too strong but this 
could be rectified by reducing the Si abundance by a further 
factor two. We suspect that this is a real effect but are hes- 
itant to be adamant on this point, given the disagreement 
between theory and observation for the CI, II lines. 

The residuals are therefore not dominated by statisti- 
cal errors but instead by systematics due to inadequacies in 
both the calibration of the data and the input physics of 
the theoretical models. Therefore the formal errors emerg- 
ing from the standard Levenberg-Marquardt algorithm will 
underestimate the errors of the atmospheric parameters. In 
addition, for white dwarfs this cool there is a well known 
correlation between effective temperature and surface grav- 
ity: a model with a given effective temperature and surface 
gravity is indistinguishable (at the level of the residuals in 
Fig. 10) in its relative flux distribution to another model 
with both a higher temperature and higher surface gravity, 
or both a lower temperature and lower surface gravity. 

In order to arrive at an error estimate, we therefore em- 
ployed the V magnitude of the white dwarf as a discrimina- 
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G140L + G230L s models 14220/8.34 13890/8.14 14550/8.56 



Table 8. White Dwarf Radial Velocities 




Orbital 
Phase 



0.767 
0.801 
0.826 
1.154 
1.187 
1.212 
1.237 
1.261 



Heliocentric 
Radial Velocity 
(km s- 1 ) 

176 
214 
253 
-71 
-80 
-78 
-76 
-80 



Orbital 

Phase 



1.286 
1.599 
1.624 
1.666 
1.690 
1.706 
1.740 



Heliocentric 
Radial Velocity 
(km s- 1 ) 

-26 
131 
115 
129 
165 
207 
234 



Filled circles-WD Asterisks-dMe Pentagons-WD+dMe 




Figure 11. Upper: sum of G140L and G230L spectra of 
EC13471-1258 on 1999 Aug 28, along with three models with 
effective temperatures and surface gravities of (13890, 8.14), 
(14220, 8.34), (14550, 8.56). Lower: as for the upper panel ex- 
cept that the wavelength range has been shifted to the optical. 
In addition, the optical spectrum from Fig. 2 is shown along with 
UBVR measurements for the white dwarf (filled circles), M dwarf 
(asterisks) and their sum (pentagons). See text for further details. 
The abscissa is in A and the ordinate is in erg s _1 cm~ 2 A -1 . 



tor. The upper panel of Fig. 11 shows a comparison between 
the sum of the G140L and G230L data and three models, 
one of which is the best fitting model (14220, 8.34) shown in 
Fig. 10. The other two models have atmospheric parameters 
of 13890 and 14550 K, with corresponding surface gravi- 
ties of 8.14 and 8.56, respectively. All three models have 
been normalized to the G140L observations at 1450 A and 
all three provide a reasonably close match over the G140L 
wavelength range (due to the correlation between tempera- 
ture and surface gravity mentioned above). Note, however, 
that the G230L data lie above all three models in the range 
1600-2150 A. We attribute this discord to STIS calibration 
uncertainties. 

The lower panel of Fig. 11 shows the same three mod- 
els as in the upper panel but extended to the optical region. 
The sum of the G230L data are shown, along with the optical 
flux distribution from Fig. 2, and the UBVRI photometry 
of the white dwarf (filled circles), M dwarf (asterisks) and 
the combination (pentagons) (Tables 4 and 7) . The photom- 
etry was converted to fluxes using the zero points in Bessell 
(1979). The best-fitting model, (14220, 8.34), passes through 
the V measurement at 5500 A with the (13890, 8.14) model 



Template 

Unshifted Sum 
Shifted Sum 
Model 



(Ki sin i) 
(km s- 1 ) 
96 ± 8 
124 ± 9 
137 ± 10 



7i,0 
(km s- 1 ) 



61 ± 10 



lying above and the (14550, 8.56) model lying below. We em- 
phasize that none of the data in Fig. 11 has been adjusted 
for consistency but simply taken from the measurements 
and associated calibrations. The ability of the photometry 
to distinguish the different models relies on the absolute cal- 
ibration of the G140L data and the V measurement. It is our 
belief that the relatively poorer fit to the V measurement of 
the (13890, 8.14) and (14550, 8.56) models represents a rea- 
sonable estimation of ±lcr for the atmospheric parameters. 
For comparison with the \ 2 technique, the models bracket- 
ing the best fit in Fig. 11 fall on the A\ 2 = 2.3 contour of the 
error ellipse (Fig. 15.6.4 of Press et al. 1992). We stress the 
correlation of the errors so that if the real effective tempera- 
ture of the white dwarf is either higher or lower than 14220 
K, the surface gravity must correspondingly be adjusted up 
or down. 

As mentioned already, the metal lines play no role in the 
effective temperature and surface gravity determination and 
as the fits were imperfect, it was deemed acceptable to judge 
the fitting by eye. A reasonable overall fit was obtained with 
abundances of 1/30 solar (except that Si was decreased by a 
further factor two); 1/10 as well as 1/100 solar models were 
clearly less satisfactory. We deduce from this a conservative 
error estimate of ±0.5 dcx. 

In summary, the modelling of the white dwarf spectrum 
has yielded the following results: 
T cff = 14220 ± 350 K 
logg = 8.34 ±0.22 
log Z = -1.5 ±0.5 



4-1-3 White Dwarf Radial Velocities 

One of the prime motivations for obtaining HST ultravi- 
olet spectroscopy was to measure the radial velocity vari- 
ations of the white dwarf arising from the binary motion. 
The fifteen G140L spectra were timed to occur as close to 
binary quadrature as possible, given the limited time avail- 
able and the constraints of the HST orbit. Radial veloci- 
ties were extracted by cross-correlating the fifteen individual 
spectra with a template (Tonry & Davis 1979) from which 
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Figure 12. White dwarf radial velocities from the HST spectra 
relative to the model template (filled circles) with best fitting 
orbital solution superimposed (solid line). 



the continuum and other low order variations has been re- 
moved. Three separate templates were tried: (i) the sum of 
all the G140L spectra; (ii) the best fitting model (middle 
curve in Fig. 10); (iii) the sum of all the G140L spectra 
but with each having been shifted to the rest frame of the 
white dwarf by the preliminary solution to the radial veloc- 
ity curve (upper curve in Fig. 10). This latter technique is 
intended to sharpen the features in the summed spectrum 
which would otherwise be slightly degraded by the orbital 
motion. For all templates and individual spectra, the con- 
tinuum was removed prior to cross- correlation. In addition, 
Lya is too broad to contain sensible velocity information 
so this line was treated as continuum and removed as well. 
The resulting spectrum's strongest feature was the quasi- 
molecular feature at 1400 A with additional velocity 
information coming from the metal absorption lines. 

The resulting velocities, for the model template only, 
are listed in Table 8 and plotted in Fig. 12. Both the motion 
of the earth with respect to the sun and the spacecraft with 
respect to the earth have been accounted for in the values 
listed in Table 8 (the latter is, of course, very small: less 
than 5 km s _1 ). 

The orbital ephemeris was used to calculate the orbital 
phase of the velocities and a sinusoidal curve with fixed or- 
bital period and orbital phase of the form: 

V r = 71 — (Ki sin i) sin(27r orb ) 

was fitted to the velocities by least squares (with only 15 
data points, the best results are obtained if the fewest num- 
ber of variables are allowed) . The best fitting curve for the 
model template velocities is plotted as the solid line in Fig. 
12. Table 8 also shows the orbital semi-amplitudes, Ki, for 
the three different templates, along with the formal errors of 
the fit. As expected, the unshifted sum template has a lower 
velocity than the other two which are consistent with each 
other within the errors of measurement. The model template 
velocities also provide a mean, or gamma, velocity for the 
fit, which we shall use in detecting the gravitational redshift 
of the white dwarf. 



Table 9. Optical Spectroscopy Observing Log (superscripts in 
column 1 refer to Table 10) 



Date 


Time 


Exp. 


Wave- 


Resol- 


No of 




Span 


Time 


length 


ution 


Spcct. 








Range 








CUT) 


(sec) 


(A) 


(A) 




1992 Apr 












02/03 


21:14-21:57 


2400 


5000-7400 


4.4 


2 




22:12-03:32 


600 


3800-5100 


2.2 


24 


05/06 


02:41-03:36 


2400 


5000-7400 


4.4 


1 


1992 May 












09/10 


18:02-01:34 


600 


5500-6700 


2.2 


34 


11 /io 
11/ 1Z 


1 7./in OO-Q c; 


yuu 


OOUU-OZUU 


o.U 


1 Q 


1993 Feb 












24/25 


22:55-02:55 


900 


3800-5100 


2.2 


16 


25/26 


22:59-02:46 


900 


3800-5100 


2.2 


15 


26/27 


22:51-02:50 


900 


3800-5100 


2.2 


14 


1993 Apr 












14/15 1 


19:14-02:02 


600 


5600-6800 


2.2 


34 


15/16 2 


21:26-21:53 


2400 


5000-7400 


4.4 


1 


15/16 


23:34-23:55 


1200 


3700-7100 


6.0 


1 


16/17 3 


18:58-01:02 


600 


5600-6800 


2.2 


19 


18/19 4 


19:00-02:27 


600 


5600-6800 


2.2 


36 


19/20 5 


19:14-00:24 


600 


5600-6800 


2.2 


26 


1993 June 












24 


10:27-14:41 


900 


5600-6800 


2.8 


14 


25 


08:17-13:30 


900 


5000-7400 


2.8 


21 


26 


08:33-13:31 


900 


3700-7100 


2.8 


18 


27 


09:06-14:18 


900 


5600-6800 


2.8 


20 



4.2 Optical Spectroscopy 

Table 9 shows an observing log for the optical spectroscopy. 
With the exception of June 1993, all the rest of the data were 
obtained with the Image Tube Spectrograph on the SAAO 
1.9-m telescope, equipped with an S20 photocathode and 
Reticon diode array detector. Much of the data was obtained 
for radial velocity measurements and great care was taken to 
ensure the stability of the wavelength scale of the spectra by 
observing Cu-Ar or He-Ne arcs every 20 min. Pixel to pixel 
sensitivity variations in the detector were calibrated using 
observations of an incandescent lamp. Spectrophotometric 
standards were also observed to achieve flux calibration but 
as these were mostly obtained through a narrow slit, the 
absolute calibration is not reliable. One of the standards 
from Hamuy et al. (1994), EG99, is only ~ 7 degrees away 
from EC13471-1258. Despite this, it is also possible that 
wavelength-dependent slit losses are present, especially at 
the extreme ends of the wavelength range covered, due to 
the fact that the spectrograph slit was aligned in the E-W 
direction and the image tube chain suffers some variability 
in sensitivity at the edges of the detector. The data were 
reduced by fiat- fielding, fitting a 5th order polynomial to the 
arc lines to define the dispersion relation, subtracting the sky 
observed in the second Reticon diode array 30 arcsec from 
the object on the sky, and flux calibrating with a relation 
derived from the standard star observations. 

Spectra were also acquired on the nights of 1993 June 
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Table 10. Observed M Dwarfs (superscripts in column 1 refer to 
Table 9) 



Name 




Spec. 




V - 1 v 


V T 

V -1 


1 LvJo 


r 1 a m 


GJ 




Type 


Vel. 

(km/s) 










514 


2 


M0.5 




0.98 


2.04 






382 


5 


Ml. 5-2 




1.00 


2.17 






393 


2,5 


M2 




1.02 


2.24 






381 


2,3,4 


M2.5 


31,33 


1.04 


2.32 






273 


1,2 


M3.5 


18 


1.17 


2.71 


0.43 


0.86 


486 


2 


M3.5 




1.17 


2.69 


0.39 


0.88 


447 


2,5 


M4 


-31 


1.29 


2.97 


0.32 


0.83 


285 


5 


M4.5 


26 


1.27 


2.95 






299 


2 


M4.5 




1.25 


2.92 






551 


1,2,4 


M5.5 




1.65 


3.65 






M dwarf 


M3.5- 




1.20 


2.75- 


0.32 


0.89 


in 13471- 


M4 






2.95 






1258 

















24-27 at Mount Stromlo & Siding Spring Observatories with 
the 2.3-m at Siding Spring, equipped with the Double Beam 
Spectrograph. Red and blue spectra were acquired simulta- 
neously with the vast majority of exposures being of 900 s 
duration but with occasional shorter exposures of 600 and 
even 450 s in good conditions. PCA detectors were used in 
each arm of the spectrograph. The red spectra covered the 
wavelength range 6420-6680 A with a resolution of ~ 2.8 A; 
the blue spectra covered the wavelength range 3900-4500 A 
with a resolution of ~ 3.0 A. 

Only one arc per night was observed. The stability of 
the wavelength scale was achieved by correcting the observed 
sky spectra so that the faint night sky emission lines were 
placed at their correct wavelengths. In the case of the blue 
spectra, the first night showed no sky lines useful for this 
purpose so these data were not used where wavelength ac- 
curacy was required. 

Whereas contributions to the spectrum in the HST data 
originate only in the white dwarf, in the optical spectrum 
both white dwarf and M dwarf contribute, in roughly equal 
amounts at 5500 A. Accordingly it is necessary to disentan- 
gle these contributions. The relatively reliable model for the 
white dwarf derived above shows a smooth flux distribution 
in the optical (Fig. 11); longward of 5500 A, its contribu- 
tion declines and Ha is its only spectral feature so we de- 
cided to concentrate on the red spectra first to measure the 
properties of the M dwarf. Our approach is as empirical as 
possible: despite substantial progress over the last decade in 
modelling all aspects of M dwarfs, there is still discord be- 
tween theory and observation (summarized in Reid & Haw- 
ley 2000). Accordingly, in addition to the spectrophotomet- 
ric standard stars, M dwarfs were also observed (Table 10), 
for the purposes of spectral type and radial velocity deter- 
mination. The superscripts in Table 10 next to the M dwarf 
star names indicate that that star was observed on the date 
in Table 9 with the same superscript. Spectral types were 
taken from Kirkpatrick, Henry & McCarthy (1991), Reid, 
Hawley & Gizis (1995) and Hawley, Gizis & Reid (1996); 
for stars with multiple spectral type estimates, the different 
estimates were identical except in the case of GJ 382 as in- 



dicated in Table 10. Radial velocities with an accuracy of 
better than 1 km s~ were taken from Nidever et al. (2002) 
except for GJ 381 which is a known binary (Delfosse et al. 
1999); Dr. Xavier Delfosse kindly supplied the radial veloc- 
ity appropriate for the system's very much brighter primary 
star (in the R band) at the time of our highest resolution 
observations: 31 km s _1 on 1992 May 09/10 and 33 during 
April 1993. 

We also include in Table 10 VRI colours and TiO and 
CaH band indices which we shall use later. The VRI data 
were taken from Koen et al. (2002) or Bessell (1990); overall 
agreement between these authors is discussed in Koen et al. 
and the specific stars observed in common were different by 
at most 0.02 mag. Photometry for one object, GJ 299, was 
taken from Table Al of Reid & Hawley (2000). The molec- 
ular band indices are defined in Table 2 of Reid, Hawley 
& Gizis (1995); the data in Table 10 are our measurements 
which were made on the night of 1993 Apr 15/16; we esti- 
mate an error of 0.04 in the indices. Comparison with the 
corresponding measurements in Reid, Hawley & Gizis (1995) 
measurements show that our estimate of CAH1 is system- 
atically larger by 0.05. 



4-2. 1 Radial Velocities of the M Dwarf 

The radial velocity semi-amplitude of the M dwarf is likely 
to be in excess of 200 km s _1 , and its absorption lines are 
also likely to be considerably broadened by rotation, so be- 
fore examining its spectrum it is necessary to measure and 
account for the blurring effect of these motions. For this 
purpose, we selected the highest resolution data among the 
SAAO red spectra, obtained on 1992 May 09/10, and 1993 
Apr 14/15, 16/17, 18/19, 19/20. 

Examination of these data showed that the radial veloc- 
ity curve of the M dwarf could be traced by the Ha emission 
line, by the Na D absorption lines or by the TiO molecular 
bands and other weaker metal lines. In the case of Ha, it is 
possible that it will not reflect the true motion of the centre 
of mass of the M dwarf because the majority of the emis- 
sion might arise on the side facing the white dwarf due to 
the heating effect of the white dwarf. We attempted to use 
the Na D lines but achieved results of marginal reliability 
due to the weakness of the lines. Instead we used the TiO 
bands and the forest of weak metal lines as the most reliable 
indicator of the motion of the centre of mass of the M dwarf. 

Anticipating the results on the rotational velocity and 
spectral type presented below, cross-correlation templates 
(Tonry & Davis 1979) were constructed over the wavelength 
range 5600-6500 A using the spectra of the M dwarfs GJ273, 
GJ381, GJ285 and GJ447 which were observed with the 
same instrumental setup (Tables 9 and 10). These spectra 
were broadened by convolution with a Doppler profile with 
V2 sin i = 125 km s~ . The overall shape of the continuum 
was removed by subtracting a 4th order polynomial: this 
procedure retained the sharp edges of the molecular bands, 
as well as the absorption lines. 

The radial velocities extracted from the five nights of 
data were corrected for the earth's motion along the line of 
sight to the program star as well as the template star (Table 
10). Orbital phases were calculated from the ephemeris in 
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M Dwarf Absorption Radial Velocities 
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Figure 13. M dwarf absorption feature radial velocities relative 
to the template GJ447 (points) on the nights of 1992 May 09/10 
and 1993 Apr 14/15, 16/17, 18/19, 19/20 (JD 2448752, 2449092, 
2449094, 2449096 and 2449097). The best fitting orbital solution 
from Table 11 is superimposed (solid line). 



Table 11. M dwarf absorption radial velocities 



Template star 


(K2 sin i) 


72,0 


GJ381 


272 ± 6 


0±4 


GJ273 


272 ± 6 


-23 ±4 


GJ447 


260 ±6 


13 ±4 


GJ285 


260 ±6 


1±4 


Mean 


266 ± 6 


-2 ±4 



Table 2 and a function of the form: 

V r = 72 + (K 2 sin i) sin(27r orb ) 

was fitted by least squares. The results for the four tem- 
plate stars are listed in Table 11 and, in the case of GJ447, 
plotted in Fig. 13 as a typical illustration of the quality of 
the results. There is a tendency in Table 11 for the earlier 
template stars to give a larger radial velocity amplitude. As 
shown below, it is difficult to pin down the spectral type 
of the M dwarf in EC13471-1258 to better than 0.5 spec- 
tral subclasses, if for no other reason than the spectral type 
varies judging from the V — I colour variation. We believe 
that the uncertainties are dominated by the choice of tem- 
plate star; we therefore adopt a conservative estimate for 
the uncertainties and adopt the same error estimate for the 
estimate of the semi-amplitude and mean or 7 velocity as 
for the individual measurements. The results are 266 ± 6 km 
s _1 and — 2 ± 4 km s _1 , respectively. 
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Figure 14. Upper: The top curve is the sum of the 2.2 A res- 
olution red spectra from 1993 Apr 16/17 and 18/19 of GJ381. 
The bottom curve is the 2.2 A resolution red spectrum of GJ273 
obtained on 1993 Apr 14/15. The middle curve is the sum of the 
2.2 A resolution red spectra of EC13471-1258, shifted into to the 
red frame of the M dwarf before summation, and with the contri- 
bution of the white dwarf subtracted. All spectra have been left 
in extinction corrected counts (to avoid amplifying the noise at 
the ends of the spectral range) and normalized to 1.0 just short- 
ward of Ho. The spectra of GJ381 and GJ273 have been shifted 
vertically upwards and downwards for clarity. Lower: as in the up- 
per curve but with the wavelength scale expanded around 6100 
A. In addition, the original spectra of GJ381 and GJ273 have 
been convolved with Doppler profiles with V sin i = 50, 100 and 
150 km s~ 1 to enable the rotational velocity of the M dwarf in 
EC13471-1258 to be estimated. 

4-2.2 Spectral type and effective temperature of the M 
dwarf 

There are a variety of approaches for determining the spec- 
tral type of the M dwarf in EC13471-1258: 

• The least squares fitting approach of Kirkpatrick, Henry 
& McCarthy (1991). We used this technique on the 4.4 A 
resolution red spectra from 1993 Apr 15/16 and determined 
that the spectral type was later than that of GJ381 (M2.5) 
and earlier than GJ299 (M4.5). It was difficult to distinguish 
among the spectral types in between. Similar conclusions can 
be drawn from the higher resolution (2.2 A) red spectroscopy 
from 1993 Apr: Fig 14a shows a comparison of the spectra 
of GJ381 (M2.5)(top), GJ273 (M3.5) (bottom) with that of 
the M dwarf in EC13471-1258 (middle). The radial velocity 
smearing and contribution of the white dwarf has been re- 



18 D. O'Donoghue et al. 



moved from the EC13471-1258 spectrum while the spectra 
of the comparison stars have been broadened with a Doppler 
profile of 125 km s _1 (see discussion below on this point). It 
is apparent that the Na D absorption is not nearly as strong 
in EC13471-1258 as in GJ381. The molecular bands in the 
latter have smaller equivalent width too. 

• The narrow band spectrophotometry technique of Reid, 
Hawley & Gizis (1995): see their fig.2 and table 2. The TI05 
index measures the band head depth of the TIO band at 
7042 A and is the best estimate of spectral type. In the case 
of EC13471-1258, this indicates a spectral type of M4.5. 
However, the index is very sensitive to the contribution of 
the white dwarf. This in turn depends on the temperature 
of the white dwarf: as the contribution of the white dwarf 
to be subtracted decreases, the index gets larger and the 
spectral type earlier. We estimate that the measured TI05 
has an uncertainty of about 0.05 due to this effect and so 
this index does not provide greater discrimination than the 
Kirkpatrick et al. (1991) technique. 

• Probably the most sensitive indicator of spectral type 
and temperature in mid M dwarfs is V — I colour. Indeed, a 
plot of this colour against spectral type for the catalog data 
presented in Reid et al. (1995) and Hawley et al. (1996) 
shows a tight linear correlation about the point (V — I, Sp- 
Type) = (2.7, 3.5) with a slope of 2.5 (i.e. stars with spectral 
type M3 (or M4) have V - I of 2.5 (or 2.9) respectively. In- 
spection of Table 10 confirms this impression. V — I colours 
for the M dwarf (i.e. with the white dwarf subtracted) were 
calculated from the CCD photometry (Table 6) yielding a 
mean value of 2.86 and a range from 2.95 to 2.75. A plot 
of the data on HJD 2449163 is shown as the lowest curve 
(crosses) in Fig. 9. From this we conclude that the spectral 
type of the M dwarf in EC13471-1258 varies between M3.5 
and M4. This conclusion is consistent with the previous two 
bullets and all other information at our disposal. It is inter- 
esting to note that the variation in V — I in Fig. 9 is much 
smaller than the modulation in V. 

From table 4.1 or fig. 4.13 of Reid & Hawley (2000), we 
derive a mean temperature for the M dwarf of 3100 ± 50 K 
with intrinsic variation of semi-amplitude 75 K about this 
mean. Whether or not this variation is simply a function of 
orbital phase, or whether it is due, for example, to slowly 
migrating spotted regions on the surface of the M dwarf, will 
require far more data than is available in this study. 



4-2.3 Rotational velocity of the M dwarf 

Fig. 14b shows the same spectra as in Fig. 14a but with 
an expanded wavelength scale and a number of different ro- 
tational models. It is clear that both rotation and spectral 
type contribute to the depths of the lines. If GJ381 is used 
as the template, a rotational velocity of 150 km s _1 is de- 
duced. However, as shown in the last section, GJ381 is too 
early. The spectral fit to GJ273 (lower spectra) in parts of 
the wavelength range covered is too poor to judge the rota- 
tional velocity. No objective technique could be devised and 
inspection was resorted to. We conclude that the rotational 
velocity of the M dwarf is V2 sin i = 125 ± 25 kms -1 . 



4.2.4 Metallicity of the M dwarf 

With considerably subsolar metal abundances in the white 
dwarf, it is of some interest to determine the composition of 
the M dwarf. As illustrated in fig. 2.18 of Reid & Hawley 
(2000), chemical composition in M dwarfs is indicated by 
the increase in strength of the CaH compared to the TiO 
molecular bands as metal abundance decreases. Quantita- 
tive measurement is described by Gizis (1997); as shown by 
his fig. la, the distinction between disk M dwarfs on the 
one hand and subdwarfs or extreme subdwarfs on the other 
(with metal abundance of 1/16 solar or less) is clearest in the 
CAH1-TI05 diagram. The measurements of these band in- 
dices rely on narrow band spectrophotometry and in order to 
reduce systematic uncertainties in our spectra of EC13471- 
1258, we measured these band strengths for three M dwarfs 
of similar spectral type (Table 10). As discussed above, for 
the M dwarf in EC13471-1258, a reasonably accurate sub- 
traction of the spectrum of the white dwarf was required 
and additional uncertainty is associated with this. 

The indices for the M dwarf in EC 13471-1258 match 
those of the similar stars (all of which are normal disk M 
dwarfs); moreover, even after correcting for the systematic 
overestimation of the CAH1 index, it clear that the indices 
for the M dwarf in EC13471-1258 lie well away from the 
metal deficient M dwarf region in fig. la of Gizis (1997). 
We conclude that the M dwarf in EC13471-1258 has nor- 
mal solar abundance and defer until later discussion of the 
considerably lower metal abundance in the white dwarf. 



4.2.5 The behaviour of the emission lines 

It was obvious even at the telescope that the profiles of the 
Balmer emission lines were comprised of more than one com- 
ponent; this constrained us from using them for radial veloc- 
ity measurements as it would not be apparent how to inter- 
pret the results. We therefore folded all the SAAO highest 
resolution red spectra into 20 phase bins using the orbital 
ephemeris and scaled each in such a manner that the contin- 
uum just outside the Ha emission line profile was normal- 
ized to unity. The results are shown as a trailed spectrum 
in the left panel of Fig. 15. The right panel shows the same 
procedure applied to the MSSSO data. Bearing in mind the 
worse resolution of the latter data, both panels are consis- 
tent with each other, indicating that the features seen are 
not transient but repeatable on a time scale of months. 

The left panel shows that the Ha emission line is in- 
deed comprised of two components, leading to a doubling 
of the line profile at orbital phase 0.25, and a single peak 
at the phases of conjunction, 0.00 and 0.50. The stronger of 
the two components is well matched to the curve labelled 
2; this curve is not a fit to this component but is instead 
the motion of the M dwarf derived from the absorption line 
radial velocity curve derived above. 

Turning to the weaker of the two components, it is ap- 
parent that this varies around the orbital cycle, peaking in 
strength at orbital phase 0.20-0.25. At orbital phase 0.75, 
although much weaker in strength than half a cycle ear- 
lier/later, it is certainly present: examination of the indi- 
vidual spectra showed a clear doubling of the emission line 
in the orbital phase interval 0.60 to 0.70, disappearance of 
the weaker component around orbital phase 0.75, a brief 
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Figure 15. Left: Trailed SAAO spectra of the Ho emission line, folded on the orbital ephemeris. Two cycles arc shown for the sake of 
continuity, the second cycle being a copy of the first. Each spectrum was scaled so that the continuum just outside the Ha profile is 
normalized to unity. The grey scale is chosen so that white corresponds to a value of 1 (i.e. the continuum) and black to 2.2 (the peak of 
the emission profiles in all spectra was 2.5). Orbital phase is shown on the ordinate axis and wavelength is shown on the abscissa. The 
curve labelled 2 corresponds to the expected motion of the M dwarf as deduced from the absorption line radial velocity curve derived 
above. The curve labelled 1 corresponds to the expected motion of the white dwarf as derived from the velocity variations in the HST 
spectra; the gamma velocity of 60 km s _1 is included in curve 1. The curve labelled 3 corresponds to the motion of the feature in the 
Dopplcr tomogram at velocity co-ordinates (75, -75) km s _1 . It can be seem as a faint feature visible between orbital phase 0.6 and 0.9 
and lying to the blue of the expected motion of the white dwarf. See text for more details. Right: As in the left panel except that the 
MSSSO data are shown. 



reappearance at orbital phase 0.80 before merging with the 
redward moving stronger component at eclipse. 

The curve labelled 1 corresponds to the expected mo- 
tion of the white dwarf as derived from the velocity varia- 
tions in the HST spectra, including its gamma velocity of 
60 km s _1 . We are not convinced that this curve is a sat- 
isfactory description of the velocity variation of the weaker 
component in the data because it co-incides with no emission 
between orbital phase 0.70 and 0.85. Omitting the gamma 
velocity (which we attribute to a gravitational redshift from 
emission near the surface of the white dwarf), makes the 
fit worse. Whatever the case, the weaker component must 
originate on the same side of the binary centre of mass as 
the white dwarf as its curvature is clearly of opposite sign 
to that of the M dwarf's motion. 

We are convinced of the reality of the feature visible 
between orbital phases 0.60 and 0.90 just slightly to the blue 
of the curve labelled 1. The curve labelled 3 is derived from 
the Doppler tomogram discussed below. Note that its motion 
does not co-incide with the arc which is strongest between 
orbital phases 0.00 and 0.25 and which the curve labelled 1 
passes through. Thus neither curves 1 nor 3 provide, on their 
own, a satisfactory representation of the weaker component. 

It is also important to note the variations in strength 
of the stronger component. At phases close to conjunction, 
there is an obvious increase in strength which is plausibly 
explained by the summation of the two components. How- 
ever, there is additional intrinsic variation in the strength 
of the stronger component: it weakens by at least a factor 
of two in the orbital phase interval 0.60-0.70 compared to 
the combined strength at orbital phase 0.50, and strength- 
ens again suddenly at orbital phase 0.75. It is clear that 



the profile cannot be explained by the simple summation 
of two optically thin emission line sources of variable radial 
velocity. Turning briefly to the blue optical spectra (Fig. 17) 
which will be discussed more fully below, the relative fluxes 
in the Balmer emission lines shows that they are optically 
thick: the Balmer decrement is flat - the line strength ratios 
Ha/H/3, H7/H/3, and US/U/3 are 1.4, 0.67 and 0.83, whereas 
Case B recombination predicts values of 2.80, 0.47 and 0.29, 
respectively (table 4 of Ferland et al. 1982). This provides a 
natural explanation for the intrinsic line strength variations 
in the Ha profile: absorption of Balmer line photons by in- 
tervening gas is easily capable of explaining the variations 
seen. Specific scenarios in which this is taking place will be 
discussed later. 

A Doppler tomogram (Marsh & Home 1988) of the Ha 
emission line in the SAAO data was calculated using the 
code of Spruit (1998). The results are shown in Fig. 16 
which includes both velocity and spatial co-ordinates (the 
former on the lower horizontal and left vertical axes, the 
latter on the upper horizontal and right vertical axes). The 
origin of the spatial co-ordinates is the white dwarf and that 
of the velocity co-ordinates the binary centre of mass (de- 
picted by a cross). The scales were adjusted, using the bi- 
nary parameters adopted from later analysis, so that unity 
in the spatial system corresponds to (Ki + K2) sin i in the 
velocity system. The centre of mass of the M dwarf and the 
white dwarf are thus at their correct position in both co- 
ordinate systems, and so are the Roche lobes. We are aware 
that this implies that the velocity of any specific feature in 
the Doppler map arises only from the binary rotation and 
no other velocity component is included (e.g. the velocity of 
any streaming motions). Indeed, in velocity co-ordinates the 
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Figure 16. Left: Dopplcr tomogram of the Ha emission line. Two 
scales are shown: velocity (lower horizontal and left vertical) and 
spatial (upper horizontal and right vertical). Both Roche lobes 
are shown as well as the binary centre of mass (cross), the white 
dwarf position (small filled circle) and the mass transfer stream 
in spatial co-ordinates as it loops around the white dwarf. The 
small circle surrounding the white dwarf is the inviscid disc radius 
(eq. 2.19 of Warner 1995). The Doppler map is scaled so that, in 
velocity co-ordinates, the M dwarf and white dwarf are at their 
correct respective positions (and by implication the Roche lobes 
are as well). The vertical line from the LI point is the mass trans- 
fer stream in velocity co-ordinates. See text for discussion of the 
significance of all these features. 

mass transfer stream, as it moves away from the LI point, 
quickly acquires ballistic velocities well in excess of the local 
binary rotation velocity and this is illustrated in Fig. 16 by 
the mass transfer stream whose velocity locus is the almost 
vertical line from the LI point. Checks on the velocity of the 
stream shows that it never intersects the Roche lobe of the 
primary star. 

With these considerations in mind, it is apparent that 
the majority of the emission arises from the vicinity of both 
Roche lobes. Further, the strongest emission component is 
centred on the Roche lobe of the M dwarf and has no bias 
towards or away from its inner hemisphere facing the white 
dwarf. This emission component corresponds to the line la- 
belled '2' in Fig. 15. 

Within the white dwarf Roche lobe, the emission is bi- 
ased towards the side facing the M dwarf and there are two 
distinct features at velocity co-ordinates (Vx, Vy) = (-25, 
-75) and (75, -75) km s _1 . The latter corresponds to the 
curve labelled '3' in Fig. 15. The former feature does not 
fit the emission in the phase interval 0.0 to 0.25 as well as 
the curve labelled T' in Fig. 15. The interpretation of these 
two distinct features in the Doppler map is not straightfor- 
ward. Neither occurs at a "special" place in the binary: e.g. 
neither feature is associated with the white dwarf. Further- 
more, it is clear that while there is no accretion disc in the 
system of the kind seen in cataclysmic variables, there may 
be a very weak mass transfer stream flowing from the LI 
point. In this case, one might expect features in the Doppler 
map corresponding to where this stream collides with it- 
self. This could lead to an inviscid (i.e. dissipationless) disc 
which would be very hard to detect (see Warner 1995 and 
references therein for a discussion of inviscid discs within 
the context of disc formation). If there is, there might be 



a feature in the Doppler map corresponding to where the 
stream collides with such a disc. The radius of an inviscid 
disc is shown in Fig. 16 by the small circle around the white 
dwarf (eq. 2.19 of Warner 1995). One of the two distinct 
velocity features mentioned occurs near to, but not at, the 
intersection of the mass transfer stream with itself or the 
radius of an inviscid disc. The other feature is near to, but 
not at the place where the ballistic stream slows down to 
low velocities at the edge of its Roche lobe before plunging 
back for a second passage around the white dwarf. 

We caution against expecting too much from the 
Doppler tomogram: the technique does not allow for in- 
trinsic variations in emission line strength due to optical 
thickness in the lines (which are clearly present); moreover, 
there is some variation in the maps produced by spectra 
from each half of the orbital cycle (not shown) . We are im- 
pressed that the emission line behaviour is not transient (it 
was seen in both the SAAO data on different occasions and 
in the MSSSO data). We believe that these results indi- 
cate the presence of streaming motions in the white dwarf 
Roche lobe implying the existence of mass transfer, albeit 
at a very low level. There is undoubtedly accretion on to 
the white dwarf from mass ejected by the M dwarf due to 
flares or its wind, i.e. not via the LI point. Such accretion 
would not cause specific and stable features in the Doppler 
map. Clearly more data with at least double the resolution 
in wavelength and binary phase are required to elucidate 
these issues further. 

In the right hand panel of Fig. 15, it is worth draw- 
ing attention to the additional emission occurring at orbital 
phase 0.35 and positioned to the red of the most red ward 
motion of the M dwarf. Although these are phase folded 
spectra, the MSSSO data have only 3-4 spectra per bin on 
average and examination of the data that contribute to this 
phase bin show a clear triple structure to the emission line. 
Subsequent spectra showed the additional peak following the 
motion of the M dwarf indicating that the origin of the flare 
is close to that of the M dwarf and therefore attributable 
to the flaring behaviour noted above. There is also evidence 
for similar behaviour (but with the additional velocity com- 
ponent occurring in the extreme blue wings of the line) in 
some of the SAAO spectra from 1993 Apr 19/20. Flaring 
in the Ha emission line has also been seen by Kawka et al. 
(2002). 

The behaviour of the blue emission lines H/3, H7 and 
H5 was examined using the SAAO data from 1993 Feb and 
1992 Apr 02/03, as well as the MSSSO data. The signal-to- 
noise in these data were generally much worse than in the 
Ha data shown in Fig. 15; despite this, the behaviour of 
these features was consistent with that seen in Ha: a two- 
component profile with obvious doubling at phase 0.25 but 
with one component much weaker half an orbital cycle later. 

The Call K line's behaviour was somewhat different in 
that it never showed doubling, was of lower velocity ampli- 
tude than the Balmer lines, and showed a strong variation 
around the orbital cycle with peak strength at orbital phase 
0.5 and essentially disappearing at eclipse. This behaviour is 
strongly suggestive that this feature is formed on the inner 
hemisphere of the M dwarf facing the white dwarf. 

Finally, the Ha emission line was used to provide a 
more accurate estimate of the rotational velocity of the M 
dwarf than the absorption line value derived above. This is 
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Figure 17. The top curve is the sum of the blue spectra from 
1992 Apr 02/03 and 1993 Feb. The bottom curve is the spectrum 
of the only M4.5eV star in the spectral atlas of Jacoby, Hunter 
k. Christian 1984. The middle smooth curve is the white dwarf 
model fit from the HST data; the remaining middle curve is the 
difference between the top and bottom curves. 



justified as it is clear from the Doppler tomography above 
that the velocity of the stronger component of the emission 
line is identical to that of the M dwarf. Fits of Doppler pro- 
files to the emission line, and allowing for the instrumental 
resolution, yielded a value of V ro t,2 sin i = 140 ± 10 km s _1 . 



4-2.6 The blue spectra 

The blue spectra from 1992 Apr 02/03 and from 1993 Feb 
were selected and summed together into 10 bins. Although 
the motion of the Balmer emission lines was obvious, no sign 
of the orbital motion of the white dwarf could be detected 
because the lines were too broad and the signal-to-noise ratio 
insufficient. Instead, the data were simply summed together 
and fluxed using the lower resolution spectrum in Fig. 2. 
The result is shown in Fig. 17 (top curve). It is clear that 
this spectrum is contaminated not only by Balmer emission 
but also by absorption features from the M dwarf, especially 
at the red end. We had no observations of M dwarfs in this 
spectral region, we could find no publicly available models 
with features to match the observations, so we made use of 
the observed spectrum of the M4.5eV star in the spectral 
atlas of Jacoby, Hunter & Christian (1984). Although its 
spectral type is 1 subclass later than that of the M dwarf 
in EC13471-1258, it was satisfactory for our purposes. A 
scaled version of this spectrum appears as the bottom curve 
in Fig. 17. The middle smooth curve shows the model for 
the white dwarf (T eff = 14420 K, log g = 8.34) along with 
the difference between the top and bottom curves. Allowing 
for subtraction uncertainties and instrumental effects, these 
two curves are a satisfactory match. 

We tried to extract an independent estimate of effec- 
tive temperature and gravity from the subtracted spectrum. 
This procedure was difficult as the emission cores had to 
be excluded and a fit to the continuum forced in several 
places. We decided that no reliable, independent estimate 
of the atmospheric parameters could be derived from this 
procedure. However, the solution was consistent within un- 
certainties with the HST value. 



Table 12. Summary of results 
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T e ff,i 


14220 


±350 


K 


log g 


8.34 


±0.22 




log Z 


-1.5 


±0.5 




Ki sin i 


137 


±10 


km s _1 


71 


61 


±10 


km s~ 1 


V ro t,i sin i 


400 


±100 


km s^ 1 


K2 sin i 


266 


±6 


km s~ 1 


72 
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km s~ 1 


V ro t,2 sin i 


140 


±10 


km s~ 1 


T e ff ,2 


3100 


±75 
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5 ANALYSIS OF THE BINARY 

The numerical results from the previous analyses are listed 
in Table 12 along with their errors. Error estimates for some 
quantities, e.g. the eclipse analysis, are not shown as they 
are so small compared to the other error estimates that they 
can safely be neglected. A subset of these parameters allow 
derivation of the fundamental properties of the binary using 
the standard equations of spectroscopic binaries to derive 
the masses and radii of the individual components: 



a 


= PorbQKl +K 2 )/2tt 


M 


= to 2 a 3 /GPib 


q 


= Ki/Ka 


Mi 


= M/(l+q) 


M 2 


= q Mi 



f(Ri,R 2 ,i) = 
Riobc,2/a = 0.49q 2/3 /(0.6q 2/3 + ln(l + q 1/3 )) 

where a is the binary separation, M is the sum of the masses 
of the components, q is the ratio of the secondary to primary 
mass and Riobe,2 is the volume-equivalent Roche lobe radius 
of the secondary (Eggleton 1983). 

Measurements in Table 12 show values for Ki sin i and 
K2 sin i, not Ki and K2. As is well known in the analysis 
of spectroscopic binaries, although a is only dependent on 
Ki sin i and K 2 sin i, M is dependent on the cube of these 
quantities. In the following analysis, therefore, uncertainties 
in the component masses and radii due to ignorance of the 
inclination and measurement errors in the radial velocity 
curves will be made explicit. The approach taken, therefore, 
was to try all integer values of i between 75° and 90° for 
the optimal values of Ki sin i and K 2 sin i in Table 12, and 
to repeat the calculations for perturbations on the velocity 
semi-amplitudes using the error estimates in Table 12. Mass- 
radius diagrams showing the results of these calculations for 
each component appear in Fig. 18. Several lines appear in 
these diagrams which will require some explanation before 
discussing their significance: 

• Consider first the lower panel (the white dwarf mass- 
radius diagram) and within this panel the line labelled with 
the numbers 266, 138. This line defines the locus of white 
dwarf mass-radius solutions for the optimal values of Ki sin i 
and K2 sin i from Table 12. Different inclination angles are 
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Figure 18. Mass-radius diagrams for the primary (lower panel) 
and secondary (upper) panel). The white dwarf mass-radius rela- 
tion from Wood (1995) is shown in the lower panel labelled with 
'Wood'. Mass-radius relations for lower main sequence stars are 
shown in the upper panel: equation 2.99 of Warner (1995) is la- 
belled 'Warner' and that derived in Clemens et al. (1998), based 
on Reid & Gizis (1997), is labelled 'Clemens'. Each sloping line 
represents a solution for the binary using radial velocity semi- 
amplitudes indicated by the pair of numbers at the end of the 
corresponding line. Each dot along each line represents a one de- 
gree change in binary inclination; one line in each panel is labelled 
with the numerical value of the inclination. The 6-sided star on 
the white dwarf mass-radius relation represents the estimate of 
log g from the model atmosphere fitting to the HST data. In the 
upper panel, the dotted line represents the limit where the mean 
density of the M dwarf would exceed the Roche lobe mean density 
associated with that orbital period, and would therefore overfill 
the lobe. The mass derived from the gravitational redshift is 0.9 
Mq; the vertical line in the lower panel is the acceptable mass 
range derived from the 1-cr error estimate on the gravitational 
redshift. See the text for details of other features and additional 
discussion. 



shown as filled circles on this line. Larger filled circles are 
shown every five degrees along with a label showing the nu- 
merical value of the inclination. This line extends from 90° 
(at the top) to 75°. 

• Lines parallel to the one just discussed are shown and 
represent solutions for the values of Ki sin i and K2 sin i 
indicated at the top of the line. In the upper panel corre- 
sponding lines are shown in the mass-radius diagram for the 



M dwarf. These loci will be referred to as "dynamic solu- 
tions" . 

• In the lower panel the line labelled 'Wood' indicates the 
mass-radius relation for white dwarfs from the finite temper- 
ature, evolutionary models of Wood (1995). On this line is 
a 6-sided star which indicates the mass and radius of the 
white dwarf deduced from the value of log g combined with 
the Wood mass-radius relation. The uncertainty on the value 
of log g is indicated by showing the mass-radius relation as 
a solid line. Beyond the 1-cr error interval, the mass-radius 
relation is shown as a dashed line. 

• We interpret the difference between the 7 velocities of 
the white dwarf and M dwarf as arising from the gravita- 
tional redshift of the white dwarf. This directly translates 
into a white dwarf mass of 0.9 Mq. The uncertainty on the 
difference allows a lower mass estimate and the allowable 
interval (at the 1-cr level) is indicated by the vertical line in 
the lower panel with the left pointing arrow attached and 
labelled 'Gravitational Redshift'. 

• As noted above, in the upper panel the dynamic solu- 
tions encompass a range of inclinations from 90° to ~ 75°. 
At lower inclinations, the secondary would exceed its Roche 
lobe and are presumably excluded solutions. The parts of 
the dynamic solutions so excluded lie above the dotted line 
in the upper panel. This explains why the range of inclina- 
tions extended from 90° to only 75°. Note that the dotted 
line is essentially a geometric constraint, being based on the 
period-mean density relation for Roche lobe filling secon- 
daries (see Warner 1995). 

• Mass-radius relations for lower main sequence stars are 
shown in the upper panel: equation 2.99 of Warner (1995) is 
labelled 'Warner' and that derived in Clemens et al. (1998), 
based on Reid & Gizis (1997), is labelled 'Clemens'. 

• The remaining three almost horizontal lines in the up- 
per panel, the leftmost of which is labelled 'V ro t,2', delineate 
a locus of M dwarf radii derived from the measurement of 
Vrot,2 sin i. The dependency on sin i is shown in the same 
manner as for the radial velocity lines which also define the 
corresponding value of M2. Three lines are shown for three 
combinations of Ki sin i and K2 sin i. 

Assuming that M dwarf radii in excess of the Roche lobe 
radius are excluded, the dynamic solutions representing the 
radial velocity semi-amplitudes and their 1 — cr uncertainties, 
as well as the eclipse analysis, restrict the masses of the com- 
ponents to 0.59-0.83 and 0.30-0.45 M for the white dwarf 
and M dwarf respectively. More precise values are available 
only if additional constraints are used. 

The additional, purely observational constraints avail- 
able are the loci of the M dwarf rotational velocity in the up- 
per panel, and the gravitational redshift in the lower panel. 
The uncertainty on the M dwarf rotational velocity is 7 
per cent or 1.5 carets on the ordinate axis of the upper 
panel. This restricts the range of allowable inclinations to 
less than 77°. The gravitational redshift constrains the al- 
lowable range for the white dwarf mass to be larger than 0.8 
M©. The dynamic solutions impose a corresponding lower 
limit on the M dwarf mass to be larger than 0.43 Mq. 

Inclusion of the Wood (1995) white dwarf mass-radius 
relation restricts the allowable inclinations to be smaller 
than 80°. 

Inclusion of the mass-radius relations for the M dwarf 
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leads to interesting conclusions: (i) the relation in Warner 
(1995), which is typical of mass-radius relations from a num- 
ber of sources, along with the assumption that the M dwarf 
radius does not overfill its Roche lobe, restricts the M dwarf 
mass to be smaller than 0.34 Mq. This is incompatible 
with the gravitational redshift constraint; (ii) Clemens et al. 
(1998) used the near infrared observations of Reid & Gizis 
(1997) to derive a mass-radius relation for M dwarfs, discov- 
ering an abrupt change in slope of the mass-radius relation 
as well as incompatibility with the mass-radius relations of 
the best theoretical models (see Reid & Hawley 2000 and 
references therein for a full discussion). The upper panel of 
Fig. 18 shows that the Clemens et al. relation, which di- 
verges from the Warner and other "standard" relations at a 
secondary mass of 0.3 Mq, are compatible with the dynam- 
ical solutions for the M dwarf in EC13471-1258. 

In arriving at a solution to the binary parameters, we 
attempted to achieve consistency with: (i) the M dwarf not 
overfilling its Roche lobe; (ii) the atmospheric parameter 
determination and Wood mass-radius relation for the white 
dwarf; (iii) the Clemens et al. mass-radius relation for the M 
dwarf (bearing in mind that there is intrinsic scatter in this 
relation); (iv) the V ro t,2 sin i determination for the M dwarf; 
(v) the gravitational redshift measurement; (vi) the dynam- 
ical solutions. In addition, we believe that the M dwarf is 
essentially filling its Roche lobe: evidence for this comes from 
the pre-eclipse dips (Fig. 8) without accompanying flare, and 
the specific features seen in the Doppler tomogram. If so, the 
solution will be on the dotted line in the upper panel of Fig. 
18. Even if this is not the case, it is clear that the M dwarf 
is very close to filling its Roche lobe. 

Accordingly, as the best "compromise" between all 
these constraints and their errors, we adopt an inclination 
of 75.5 ± 2.0°, a white dwarf mass and radius of 0.78 ± 0.04 
M Q and 0.011 ± 0.01 Rq, and an M dwarf mass and ra- 
dius of 0.43 ± 0.04 Mq and 0.42 ± 0.02 R . We regard the 
error estimates as conservative. It is an endorsement of stel- 
lar evolution observation and theory that all the relevant 
constraints can be satisfied simultaneously. 

5.1 Comparison with the results of Kawka et al. 

Recently, Kawka et al. (2002) have reported observations 
and analysis of EC13471-1258 along with two other DA+dM 
systems. In general, their observational results broadly agree 
with those presented here. In the following parentheses, the 
results of Kawka et al. appear first, those presented in this 
paper second: the values of white dwarf log g (8.26, 8.34), 
T eff (14080 K, 14220 K),andVi, rot sin i (140 km s"\ 140 km 
s _1 ) are all consistent within the errors (though we are scep- 
tical of their small error of 0.05 on log g, given the analysis 
of the error on this quantity shown above). There is, how- 
ever, discord in the M dwarf radial velocity semi-amplitude 
(241 ±8.1, 266 ±6). We are inclined to prefer our value as it 
was derived from the absorption features of the M dwarf and 
is also consistent with the motion of the stronger component 
of the emission lines. 

The analysis of Kawka et al. is somewhat flawed in that 
they use the V sin i value to derive the absolute radius of 
the M dwarf without dependence on sin i. They then use 
the duration of the white dwarf eclipse along with the ab- 
solute radius of the M dwarf to derive the inclination (i.e. 



neglecting that this relationship also includes a, the orbital 
separation). Fortuitously, their inclination of 73.5° is reason- 
able. 

However, this inclination and, especially, the low value 
for the M dwarf radial velocity semi-amplitude yields a large 
value of 0.58 M Q for the M dwarf mass. This, along with the 
M dwarf radius of 0.42 R Q can be excluded for two reasons: 
(i) the radial velocity semi-amplitude for the white dwarf 
implied is about 185 km s _1 (not 165 km s~ x as claimed 
by Kawka et al.). This is excluded by the results in Table 
8 and Fig. 12. Indeed, Kawka et al. derived a white dwarf 
radial velocity semi-amplitude of 155 km s _1 but with large 
error (35 km s _1 ), presumably because it seems from their 
fig. 11 that they used only one G140L HST spectrum from 
each orbit; (ii) the M dwarf mass and radius fall far from 
the locus of these quantities, either theoretical or empirical, 
derived from a large sample of M dwarfs (fig. 3 of Clemens 
et al. 1998; Reid & Gizis 1997). In addition, the spectral 
type implied by the M dwarf mass is M0 or earlier (see Reid 
& Hawley 2000) and discordant with the spectral type of 
M3.5-4 derived in this paper, or M2-M4 listed by Kawka et 
al. themselves. We believe that the M dwarf mass of Kawka 
et al. should be disregarded. 



6 DISCUSSION 

6.1 Evolutionary status of the binary 

The canonical theory for the formation of cataclysmic vari- 
ables involves common envelope evolution enforcing dra- 
matic orbital shrinkage with a hot white dwarf and an M 
dwarf emerging from the envelope of the AGB predecessor 
of the white dwarf. Following this, further orbital shrink- 
age takes place via angular momentum loss by the M dwarf 
through a weak wind with eventual contact. Further evolu- 
tion of the semi-detached cataclysmic variable requires on- 
going angular momentum loss to ensure ongoing Roche lobe 
contact by the M dwarf. These scenarios are discussed in 
Warner (1995) and references therein. 

There is also an accumulation of evidence that the mass 
transfer rate in cataclysmic variables varies widely and per- 
haps cyclically on short time scales. Mass transfer may even 
shut off and the system "hibernates", a term originated to 
explain the very existence of classical novae: in order to get a 
thermonuclear detonation with the vigour of those observed, 
the mass accretion rate was required to to be very low im- 
mediately before the nova explosion (Shara et al. 1986). The 
low mass accretion rate is now no longer believed to be so 
important (Prialnik & Kovetz 1995). More recent and direct 
evidence for variations in mass transfer rate is provided by 
the very wide range of mass transfer rates at a given orbital 
period evident in diagrams such as fig. 9.8 of Warner (1995). 
One hypothesis to explain such mass transfer rate variations 
is irradiation of the M dwarf by the white dwarf. 

In terms of the above, therefore, EC13471-1258 is either 
a post common envelope binary just at the phase where its 
M dwarf comes into contact with its Roche lobe for the 
first time, or it is a system that has been in contact and is 
currently hibernating. We believe that two pieces of evidence 
favour the latter interpretation: 

• the rotational velocity of the white dwarf: with V ro t,i = 
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400 ± 100 km s _1 , the star rotates in 120 s. Such a rapid 
rotation rate implies that there has been an earlier phase 
of spin-up by mass transfer. Indeed, inspection of table 4 
of Szkody et al. (2002a) shows that this rotation velocity is 
typical of cataclysmic variables. 

• the behaviour of the Ha emission line profile: although 
the interpretation of the Doppler tomogram is not clear, 
what is clear is that there are two or more components in 
the emission line and intrinsic variability as well, even in the 
component associated with the M dwarf. As the Balmer lines 
are optically thick, it is clear that the intrinsic variability 
can be understood by absorption by gas in the Roche lobe 
of the primary. The fact that the variation is a function of 
binary phase indicates that the absorption sites are fixed in 
the rotating reference frame of the binary and hence that 
there is a very weak mass transfer stream. 

We therefore conclude that EC 13471-1258 is a hibernating 
cataclysmic variable. 

Interestingly, the white dwarf temperature in EC13471- 
1258 is cool compared to white dwarfs in cataclysmic vari- 
ables of similar orbital period. Indeed the temperature of 
this white dwarf is more typical of systems below the period 
gap (fig. 3 of Szkody et al. 2002a) . It is close to the temper- 
ature of the pulsating white dwarf in GW Lib (Szkody et al. 
2002b), but shows no oscillations. 

The white dwarf radius and model atmosphere imply a 
distance to the system of 48 ± 5 pc, well within the reach 
of ground-based parallax programs. An accurate parallax 
would enable a more stringent test of the white dwarf mass- 
radius relationship than presented in this paper. An absolute 
V magnitude of the white dwarf was derived: 11.74. 

6.2 The M dwarf 

The analysis of the M dwarf presented in this paper shows 
that its properties are compatible with field stars of sim- 
ilar spectral type. This M dwarf is of particular relevance 
to studies of M dwarfs of the kind conducted by Reid & 
Gizis (1997) and Clemens et al. (1998) as its properties (V-I 
colour, mass and radius) place it in the segment of the M 
dwarf sequence where an abrupt change in slope of many 
stellar parameters occurs. This abrupt change in slope is 
discussed in Clemens et al. (1998) and Reid & Gizis (1997) 
to which the reader should turn for further information. 

The M dwarf exhibits the properties expected of a sec- 
ondary star in a cataclysmic variable: it is extremely ac- 
tive, presumably because of its very short rotation period. 
The erratic behaviour of the binary period is also typical 
of cataclysmic variables and presumably arises from mag- 
netic cycling in the M dwarf (e.g. Warner 1995 and refer- 
ences therein). If flares of the kind described in this paper 
occurred on the secondaries of cataclysmic variables with 
higher mass transfer rate, they would be indistinguishable 
from the flickering in the system. 

Compatibility with the Clemens et al. mass-radius rela- 
tion implies that the M dwarf is 250 K hotter than the 3100 
K derived above: see fig. 4.13 of Reid & Hawley (2000). As 
they discuss, the temperature scale for M dwarfs is uncer- 
tain by this amount. Resolution of this discrepancy must 
await further developments in M dwarf model atmosphere 
research. 



A puzzling feature of the photometry of the M dwarf is 
the fact that the ellipsoidal variation is not constant and is 
often asymmetrical with respect to the conjunction of the 
two stars. When this happens, its minimum occurs slightly 
before orbital phase 0. If this behaviour could be attributed 
to star spots on the M dwarf, it would not be puzzling. How- 
ever, spots on other chromospherically active stars migrate. 
In this case, 10 years of eclipse monitoring has shown that 
the asymmetry is either absent or always in the same sense. 

Using the distance derived from the last subsection, the 
absolute V magnitude of the M dwarf is 11.82. 

6.3 Chemical composition of the binary 
components 

The chemical compositions of the binary components are dif- 
ferent: that of the M dwarf is normal solar, that of the white 
dwarf is 1/30 solar. As has been well known for decades, the 
chemical composition of the photospheres of typical single 
DA white dwarfs is extremely metal deficient (< 10~ 6 ) due 
to gravitational settling of heavy elements in the very high 
surface gravity of the white dwarf. The presence of met- 
als in the white dwarf in EC13471-1258, therefore, is due 
to mass accretion from the nearby M dwarf (whether this 
mass transfer is from the wind of the M dwarf or Roche lobe 
overflow is not relevant to this discussion). Metals are also 
seen in the white dwarf primaries of cataclysmic variables 
where sub-solar metal abundances are common (see table 4 
of Szkody et al. 2002a). The abundances nevertheless vastly 
exceed the abundances in single white dwarfs. 

If the chemical composition difference is in steady state, 
and the mass accretion rate on to the white dwarf can be 
estimated, the diffusion coefficients for the chemical species 
visible in the HST spectra can be deduced. This is beyond 
the scope of this paper but is clearly an interesting avenue 
for future work as this system provides a natural diffusion 
laboratory. 
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